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The interplay between spin and charge transport in electrically and magnetically inhomogeneous 
£SJ ■ semiconductor systems is investigated theoretically. In particular, the theory of spin-polarized bipo- 

lar transport in magnetic p-n junctions is formulated, generalizing the classic Shockley model. The 
theory assumes that in the depletion layer the nonequilibrium chemical potentials of spin up and 
£N| ■ spin down carriers are constant and carrier recombination and spin relaxation are inhibited. Under 

the general conditions of an applied bias and externally injected (source) spin, the model formulates 
analytically carrier and spin transport in magnetic p-n junctions at low bias. The evaluation of 

■ the carrier and spin densities at the depletion layer establishes the necessary boundary conditions 
' for solving the diffusive transport equations in the bulk regions separately, thus greatly simplifying 

the problem. The carrier and spin density and current profiles in the bulk regions are calculated 
and the I-V characteristics of the junction are obtained. It is demonstrated that spin injection 
through the depletion layer of a magnetic p-n junction is not possible unless nonequilibrium spin 
I ' accumulates in the bulk regions-either by external spin injection or by the application of a large 

, bias. Implications of the theory for majority spin injection across the depletion layer, minority spin 

■ pumping and spin amplification, giant magnetoresistance, spin-voltaic effect, biasing electrode spin 
' injection, and magnetic drift in the bulk regions are discussed in details, and illustrated using the 

£C) , example of a GaAs based magnetic p-n junction. 

m ■ 
o ■ 

CN ! I. INTRODUCTION 

O . 

Active control of spin in semiconductors§ is projected to lead to significant technological advances, most importantly 
in digital information storage and processing, magnetic recording and sensing, and quantum computing.1313 Using 
semiconductors for spintronic applications-where spin, in addition to charge, is manipulated to influence electronic 
properties-has several advantages. First, integration of spintronics with traditional semiconductor technology calls 
for employing semiconductors (rather than metals) as media for spin control. Second, semiconductors are versatile 
materials, not only for their electrical properties, but also for their spin/magnetic characteristics. Doping control of 
electrical and magnetic properties, optical spin orientation and detection, bipolar (electron and hole) transport, and 
interface properties (charge and spin accumulation and depletion) leading to device concepts from p-n junction diodes 
to field-effect transistors, are among the great advantages of semiconductors over other candidates for spintronic 
materials. By allowing for the active control and manipulation of carrier spin and charge by electric and magnetic 
fields as well as by light, semiconductor spintronics creates the potential for an integrated magneto-optoelectronics 
technology. 

A generic semiconductor spintronics scheme involves three steps: injection of nonequilibrium spin into a semicon- 
ductor, spin storage, manipulation, and transfer, and spin detection. Spin injection was historically first accomplished 
optically, by illuminating a semiconductor with circularly polarized light-the so called spin orientation.Q Electrical 
spin injection (that is, spin injection from a magnetic electrode, often called simply spin injection) into semicon- 
ductors, while predicted theoretically already in the 70s,El has been-denapastrated only recentlyjj and realized as an 
injection frnm.a magnetic semiconductor ,uaa a ferromagnetic metalJ130t3 and a ferromagnetic metal/tunnel barrier 
contactBtSEM3 

Once injected, nonequilibrium spin survives for a reasonably long time when compared to typical relaxation times 
of momentum and erLetgy of the injected carriers. Room temperature spin relaxation times in semiconductors are 
typically nanosecondstata (compared to sub picosecond time scales for momentum and energy relaxation). Similar in 
magnitude are only carrier (electron and hole) recombination times, which are usually between micro- to nanoseconds. 
If not in the ballistic regime, transport of spin in a semiconductor can be characterized as carrier recombination and 
spin relaxation limited drift and diffusion. Spin typically diffuses over micron distances from the point of injection, 
sufficient for microelectronics applications. Application of large electric fields can further drag the injected spin over 
several microns at low temperatures, as in intrinsic GaAs,E3 and even up to 100 /im in n-doped GaAs.EHI (As far as 
the spin diffusion length is concerned, metals have an advantage: because of the large Fermi velocity, spin diffusion 
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lengths in metals can be as large as centimeters.) Important for device applications are studies of spin transport 
in inhomogeneous semiconductors. It has .-already been shown, for example, that spin phase can be preserved in 
transport across heterostructure interfaces,Eil that electron spin can be controlled by bias in semipagnetic resonant 
tunneling diodespj and that spin can tunnel through the transition region of tunnel diodes.E3o The final step of 
a generic spintronics scheme is spin detection. Traditionally spin in semiconductors has been detected optically 
by observing circular polarization of the recombination light u Efforts to electrically detect nonequilibri um , s pin in 
semiconductors rely, on spin-charge coupling, realized either as spin-dependent Schottky barrier transportEffij or as 
magnetoresistanceEj and galvano-voltaicEa effects. nn 

After the discovery of ferromagnetism in III-V semiconductor compounds ,cj'EJ the great-push for semiconductor 
spintronics came with the fabrication of (Ga,Mn)As which is ferromagnetic above 100 KBEj Ferromagnetic semi- 
conductors can not only serve to inject and detect spin in all-semiconductor spintronic devices, but can also form a 
basis for nonvolatile memory, opening prospects of integrated, single-chip memory and logic amplications (feasibilitiz-pf 
such prospects has been demonstrated by controlling semiconductor ferromagnetism opticalljcJ'O and electricallyEj) . 
There is a steady increase in the, number of available ferromagnetic semiconductors, includingj-a, first ecoup-IV com=. 
pound GeMnB (In,pa,Mn)As,EZI reported room temperature ferromagnets Mn-doped CdGeP,E3 GaN,E3 and GaP,c£l 
and Co-doped Ti0 2 B 

Closely following the experimenteJ [ p | nocMss. r ttia4or theoretical efforts have been dedicated to understanding electrical 



spin injection into seinicondjietoi: 



semiconductorsMBBHBBBB 



and investigating fundamental issues of spin-polarized transport in 



Another direction for fundamental spjntronics theory has been predicting and 
analyzing various spintronics device architectures for possible technological^ applications. The common goal of these 
studies is devising spin valves and structures (typically including one or several magnetic., layers) with maximized 
magnetoresistance. To this end various spin field-effect transistors have been proposedjESOEJ where the source and 
drain are ferromagnetic electrodes serving to inject and detect spin which is transported in a (typically) nonmagnetic 
channel. Spin and charge transport in the channel are controlled by ga te , biaS | th ffi uc^ | tb|C R ashba effectl 
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detectorsJZj resonant tunnelina, ipolar magnetic diodesjlj quantum-interference mesoscopic schemealjO 



proposed ^intronics device schemes include heterostructure spin ^tersSESSO'O'EaEalilJO'O and spin polar,i 
and various spin emf sources 



We have recently proposed two spintronics device schemes that take advantage of bipolwp,U^Ge.tTon and hole) nature 
of transport .in inhomogeneously doped semiconductors: a spin-polarized p-n junctionEjl£3S and a magnetic p-n 
junctionO'Lj A spin-polarized p-n junction is a p-n junction with a source spin injected externally into one or both 
regions (p and n). The source spin can be injected either optically or electrically. We have demonstrated that 
nonequilibrium spin can be injected (transfered) very effectively across the depletion layer (space-charge region), from 
both regions: by the majority carriers into the respective minority region, and, vice versa, by the minority carriers 
into the respective majority region. Spin injection (throughout the paper "spin injection" will mean spin injection 
through the depletion layer, while externally injected spin will be referred to as "source" spin) by the minority carriers 
leads to spin accumulatioiijiiji the majority region, with an effect of amplifying the spin and significantly extending the 
spin diffusion/drift length.E3 We have also shown that nonequilibrittm spin can be stored and manipulated in a spin- 
polarized p-n junction by external bias-a spin capacitance effectS Furthermore, a spin-polarized p-n junction can 
generate spin-polarized currents as a spin solar cellH when illuminated by circularly polarized light, a spin-polarized 
current flows in a p-n junction. 

Magnetic p-n junctionsEj offer even more functionality by coupling equilibrium magnetism and nonequilibrium 
spin. A magnetic p-n junction is formed by doping a p-n junction with magnetic impurities, differently in the p and 
n regions. Magnetic impurities induce large g factors of— the mobile carriers, thus the application of a magnetic field 
results in a significant spin splitting of the carrier bandsHil If the doping is so large as to induce a ferromagnetic order, 
the splitting appears also without magnetic field. The important question, of whether spin can be injected by the 
majority carriers from the magnetic majority region into the nonmagnetic minority one, was answered negative. We 
have demonstrated that only if nonequilibrium spin is generated first in the majority region, it can subsequently be 
injected through the depletion layer. Spin can be also injected through the depletion layer at large biases, since then, 
without any external spin source, nonequilibrium spin is generated by the strong electric field in the bulk regions. We 
have also shown that magnetoresistance of a magnetic p-n junction increases exponentially with increasing magnetic 
field (that is, spin band splitting) at large fields. Magnetic p-n junctions exhibit even giant magnetoresistance, when 
source spin is injected into the majority region. We have also predicted a spin- voltaic effect (the phenomenon related 
to the Silsbee- Johnson spin-charge couplingEJEJ) where charge current (or voltage in an open circuit) arises solely 
due to a nonequilibrium spin maintained in proximity to the magnetic region. Magnetic p-n junctions can also serve 
as spin valves, since the direction of the zero-bias current can be reversed by reversing either the polarization of the 
source spin or the direction of the applied magnetic field. ___ 

We have studied spin-polarized and magnetic p-n junctions mainly numerically,E3£jO by solving a self-consistent 
set of recombination- relaxation and drift-diffusion equations, and Poisson's equation. We have obtained solutions 
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for the carrier and spin densities and currents for small and large biases, and different values of magnetic fields 
and the externally injected spin polarization. Numerical solution is indispensable at large biases (large injection), 
where analytical methods are not available. Large bias solutions describe carrier and spin transport as both drift and 
diffusion, since drift currents due to electric fields are significant even outside of the depletion layer. The low injection 



regime is tractable analytically. In Rcf. 54 we have introduced a heuristic analytical model which accounts well for 
the numerical findings, and explains all the important qualitative features of magnetic p-n junctions. In fact, our 
numerical solutions show that the most interesting and potentially important properties of magnetic junctions are at 
small biases; Jarge biases may still be useful for injecting spin across the depletion layer, or extracting spin from the 
bulk regions ,E-3 as described in Sec. [VE.) 

In this paper we formulate a general model of magnetic p-n junctions (the model includes spin-polarized p-n nonmafe 
netic junctions as a particular case), following the classic formulation of Shockley of ordinary bipolar j unctions e3e3 
The model describes magnetic p-n junctions at small biases (low injection), with arbitrary external (source) spin 
injection and band spin splitting (magnetic field), within the limits of nondegenerate carrier statistics. The paper 
has the dual role of describing the fundamental properties of spin-polarized transport in inhomogeneous magnetic 
semiconductors, while presenting a model simulation, based on the recombination and relaxation limited bipolar drift 
and diffusion, of novel microelectonics spintronic devices. If semiconductor spintronics is to become a reality, then 
detailed transport analyses of the type presented in this paper are essential. The fully analytic nature of our theory 
makes our model simulation particularly useful. 

The paper is organized as follows. Section II introduces the model and formulates its assumptions and approx- 
imations. Section III describes the spatial profiles of the carrier and spin densities in the bulk regions, gives the 
boundary conditions for the densities, and discusses the I-V characteristics of magnetic p-n junctions. In Sec. IV we 
apply our theory to several cases of interest: spin injection-through the depletion layer-by the majority carriers, spin 
pumping and spin amplification by the minority carriers, source spin injection by the biasing electrode, spin injection 
and extraction at large biases, and magnetic drift effects in the carrier and spin transport. Finally, we summarize our 
findings in Sec. [v[ where we also outline strategies for applying our theory to more realistic materials structures and 
more complex spintronic devices based on magnetically inhomogeneous semiconductors. 



II. MODEL 



The basis for our model is a semiconductor p-n junction in which carrier bands are inhamogeneously spin split: 
there is a finite equilibrium spin polarization of the carriers, different in the p and n regions.E3 Large (comparable to 
the thermal energy) spin splitting of carrier bands can arise as a result of doping with magnetic impurities (which 
may, but need not, contributq-to the carrier densities). Magnetic impurities can significantly increase the carrier g 
factors (usually up to g 200E-3), so that the application of a magnetic field B induces large spin Zeeman splitting, 
2C = 9^bB, of the bands (/is is the Bohr magneton). Inhomogeneous spin splitting can be realized either by 
inhomogeneous magnetic doping in a homogeneous magnetic field, or by a homogeneous magnetic doping in an 
inhomogeneous magnetic field, or both. Our model applies equally well to ferromagnetic p-n junctions, where bands 
are spin split even at zero magnetic field. To keep the discussion transparent and to avoid complex notation, we 
consider only the conduction band to be spin split (that is, only electrons to be spin polarized), keeping holes 
unpolarizcd. This simplification does not affect our conclusions, as electron and hole transports are fully separated 
in our model. (Spin polarization of holes is treated in Appendices [a| and |b[) The layout of a magnetic p-n junction 
is shown in Fig. Ill The semiconductor is p-doped with N a acceptors (per unit volume) along the x axis from — w p 
to 0, and n-doped with Nd donors from to w n . The depletion layer forms at (—d p ,d n ). We are not concerned 
with the transition region itself- we simply assume that it is steep enough (in fact, that it changes over a region 
smaller than the Debye screening length) to support space charge, and that all the spin splitting changes occur only 
within the transition region, being constant in both p an d n re gions (the case of magnetic drift where the splitting is 
inhomogeneous also in the bulk regions is treated in Sec. IV F). 

We denote the electron density as n — n(x) and the hole density as p = p{x). The corresponding equilibrium 
values are hq and p$, and the deviations from the equilibrium values are 5n = n — uq and Sp — p — po. Electron spin 
density s = s(x) (in equilibrium sq and deviation 5s = s — sq) is a difference between the densities of spin up and spin 
down electrons: s = n-j- — nj,. As a measure of spin polarization we use the spin polarization of the carrier density 
(not current): a — s/n (in equilibrium a$ and deviation 5a — a — ao). The equilibrium properties of magnetic p-n 
junctions are discussed in Appendix |A|, where uq, po, sq, and the built-in potential Vb are calculated. The transport 
parameters of the carriers are diffusivities D nn and D np of electrons in the n and p regions, electron lifetime r np in 
the p region, and electron spin lifetime T\ p and T\ n in the p and n regions. The unpolarized holes are characterized 
by D pn and r p „, diffusivity and lifetime in the n region. Throughout the paper, unless explicitly specified otherwise, 
a single subscript denotes the region or boundary (p, n, L, or R), while a double subscript denotes first the carrier 
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type or spin (p, n, or s) and then the region or the boundary (for example, r pn is the lifetime of holes in the n 
region). Terms "majority" ("minority") will refer to electrons in the n (p) region, and similarly for holes, and not 
to the more (less) populated spin states, as is usual in the physics of magnetotransport. Similarly, the term bipolar 
bears no relation to spin, describing only the transport carried by both electrons and holes. Finally, terms "bulk" 
and, equivalently, "neutral" will denote the regions outside the depletion layer, where, at low biases, charge neutrality 
is maintained. The notation is summarized in Table |]. 

The junction is driven off equilibrium by applying bias and injecting source spin. We place contact electrodes at 
x = —w p and x = w n . We keep the left electrode general, capable of injecting electrons, Sn p = Sn(—w p ) ^ 0, and 
spin, 5s p = Ss(—w p ) 7^ 0. This boundary condition covers magnetic diodes (Ohmic contact, Sn p = 0), and magnetic 
solar cells and junction transistors (5n p ^ 0). The right electrode is assumed to be Ohmic, Sp n = 5p(w n ) = 0, but 
able to inject spin, Ss n = Ss(w n ) 7^ 0. The majority carriers in both regions are assumed constant: p = N a in the p 
side and n — Nd in the n side. The source spin injection, here considered to take place geometrically at the contacts, 
can be realized either by the contact electrodes themselves (if the electrodes are magnetic), by optical orientation 
close to the contact, or by electrical spin injection from a third electrode (say, transverse to the junction current). 
Different cases mean different boundary conditions for spin. For no w we assume a third terminal injection so that 
Ss p and Ss n are free parameters of the model; we will later, in Sec. IV D consider the case of the contact (biasing) 
electrode source spin injection, where ds n will depend on the charge current in the junction. 

To reduce the initial drift and diffusion transport problem to a simple diffusion problem in the neutral regions we 
need to know the boundary conditions for the bulk regions at the depletion layer, that is, the carrier and spin densities 
riL = n(—d p ), sl = s(—d p ) at the left (L) and ur = n(d n ), sr = s(d n ) at the right (R) boundary of the depletion 
layer. We will calculate these boundary densities in the subsequent sections. 

We use several approximations to solve our model. First, we consider only low biases, meaning that the applied 
forward voltage V is smaller than the built-in field Vf,, which is typically about 1 eV. At small biases the densities of 
the minority carriers are much smaller than the densities of the corresponding majority carriers (the small injection 
limit), the electric field is confined to the depletion layer, and the bulk regions can be considered neutral. We next 
assume that the temperature is large enough for the donors and acceptors to be fully ionized, so that n — and 
p = N a in the respective majority regions, and the carriers obey the nondegenerate Boltzmann statistics (limiting 
doping densities to about f0 18 /cm 3 for typical semiconductors at room temperature). Finally, we consider only 
moderate spin splittings (much smaller than the built-in field), perhaps no greater than 5/cbT, since greater splittings 
can severely affect the band structure, and reduce the effective band gap. 

We have also made simplifying assumptions as to the band structure of the magnetic semiconductor. First, we 
neglect possible orbital degeneracy of the bands, and treat the spin states as spin doublets. We also neglect the 
effects of magnetic doping on the band structure (that is, changes in rii, additional band offsets, etc.) and that of the 
carrier density on the band spin splitting. The latter can be important in ferromagnetic semiconductors. However, 
since it is the minority carriers which determine the transport across p-n junctions, it is unlikely that a variation in 
the carrier density would appreciably affect our conclusions. We also assume that momentum and energy relaxation 
proceeds much faster than carrier recombination and spin relaxation, so that nonequilibrium, spin-dependent chemical 
potentials describe well the junction under an applied bias and with a source spin. Finally, we do not consider orbital 
effects due to the applied magnetic field, although these can be included in our theory simply by allowing for a 
magnetic dependence of diffusivities. 



III. CARRIER AND SPIN TRANSPORT IN THE NEUTRAL REGIONS 



The transport of carriers and spin in magnetic p-n junctions can be realistically described as drift and diffusion, 
limited by carrier recombination and spin relaxation. The transport equations were introduced in Ref. ^4], and have 
been solved numerically for a few important cases in Refs. p2|]53| , |54| . Denoting the carrier (here electron) and spin 
currents as J n and J s , the drift-diffusion equations are 

J n = D n (n<f>' t + sC -n'), (1) 
J s = D n (sti+n(' -s'). (2) 

Here <pt is the total local electrostatic potential, comprising both the built-in potential 4>b and applied bias V (the 
electric field is E = —(f>' t ), and the magnetic drift is proportional to the spatial changes in the band spin splitting, 
£' (see Fig. []]). Throughout this paper we express the potentials and the energies in the units of ksT/q and fc^T, 
respectively (fcs is the Boltzmann constant, T is the temperature, and q is the proton charge). In a steady state 
carrier recombination and spin relaxation processes can be expressed through the continuity equations for electrons 
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FIG. 1: Schematics of a magnetic p-n junction. The junction is p-doped from —W p to and n-doped from to w„. The depletion 
layer (space-charge region) forms at — d p < x < d n . The upper figure depicts an inhomogeneously spin-split conduction band 
and a valence band without the spin splitting. The conduction band spin splitting in the n region is 2£ nn , in the p region it is 
2Cnp- The greater the £ is, the more is the lower band (here called the spin up sub-band) populated. The intrinsic built-in field 
across the depletion layer is Vb- For electrons the built-in field becomes explicitly spin dependent: Vbf — Vt + Cnn — C«p an< i 
Vbi = Vb — Cnn + Cnp- The lower figure depicts regions with distinct transport characteristics: CDR are the (minority) carrier 
diffusion regions and SDR are the spin (here only electron) diffusion regions. The characteristic sizes of the regions are given 
by the corresponding diffusion lengths, as indicated. The unshaded areas are the homogeneous regions, where carrier and spin 
densities assume their equilibrium values. The known (input) densities of the model are n p , s p at — w p , and Nd, s„ at w n , while 
the densities at the depletion layer, til and si on the left side and tir = Nd, sr on the right side, are calculated in the text. 
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Nd — donor density 

N a = acceptor density 

n = electron density (equilibrium no) 

p — hole density (po) 

s = spin density (so) 

a — s/n, spin polarization (qo) 

8s — s — aon, the effective nonequilibrium spin density 

J„ — electron particle current 

j„ = —qJn, electron charge current 

J p — hole particle current 

j p = qJp, hole charge current 

J s = spin current 

n p — electron density at x — —w p (no p ) 

til — electron density at x = —dp (not = nop) 

s p = spin density at x = —w p (so P ) 

s n = spin density at x = w n (son) 

sl = spin density at x = —d p (sol = sop) 

s R = spin density at x = d„ (s 0R = s „) 

Wp — Wp — dp, effective width of the p region 

w n = w n — d n , effective width of the n region 

Dnn = electron diffusivity in the n region 

D np = electron diffusivity in the p region 

T np = lifetime of electrons in the p region 

T P n = lifetime of holes in the n region 

L np — J -Dnp i"np, electron diffusion length in the p region 

Lpn = y/ Dp n Tp n , hole diffusion length in the n region 

Tip — intrinsic spin lifetime in the p region 

Lip — y/ D n pTip intrinsic spin decay length in the p region 

1/Tsp = l/r n p + 1/Tip, spin decay rate in the p region 

L sp = y/ D n pT a p, spin diffusion length in the p region 

Ti n = spin lifetime in the n region 

L sn = V DnnTin, spin diffusion length in the n region 

Vb = built-in potential 

V = applied bias 



TABLE I: Summary of the notation used in the text and in Table |l| All the spin parameters (spin density, spin lifetime, etc.) 
relate to electrons. In the brackets are the equilibrium densities. 
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(3) 
(4) 



where r is the electron-hole recombination rate and s = s — aon, expressing the fact that intrinsic-spin relaxation 
processes (spin-flip scattering, say, by phonons or impurities) conserve the local carrier density.acJ Electron-hole 
recombination also degrades spin, the fact reflected in the first term of Eq. (|J) . Equations (Q) , (0) , (g) , and (|j) , together 
with Poisson's equation </>" = —p^q/eksT), where p is the local charge density and e is the semiconductor's dielectric 
constant, fully describe the steady-state carrier and spin transport in inhomogeneous magnetic semiconductors O In 
the rest of the paper (except for Sec. IV F| ), magnetic drift force will play no explicit role, since we assume that the 
magnetic doping is uniform in the bulk regions. The inhomogcncity in the spin splitting, which is confined to the 
depletion region, will appear only through the boundary conditions. 

At low biases, the case most important for device applications, the problem of the carrier .and, spin transport in 
magnetic p-n junctions reduces to the problem of carrier and spin diffusion in the neutral regions.EZIH This observation, 
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to be useful, needs to be furnished with the boundary .conditions for the carrier and spin densities at the depletion 
layer boundary (ul, sl, tir, and sr). Shockley's modelEll evaluates the carrier densities in unpolarized p-n junctions 
from the assumption that a quasiequilibrium is maintained in the depletion layer even at applied (low) biases. This 
assumption alone is insufficient to obtain both carrier and spin densities in a spin-polarized magnetic junction. We 
use, in addition, the continuity of the spin current in the depletion layer to calculate the densities. A simple version 
of this model was introduced in Ref. |34| where it was assumed that (1) at a forward bias and with a source spin 
injected into the majority region (Ss n ^ 0) the spin current at the depletion layer, J s r, vanishes, and (2) at a reverse 
bias, and with spin injected into the minority region (Ss p ^ 0), all the spin entering the depletion region is swept by 
the large built-in field to the majority side. Assumption (1) explains spin injection of nonequilibrium spin through 
the depletion layer, while (2) explains spin pumping by the minority carriers. Both assumptions will follow as special 
cases of the spin current continuity, in our moidel. 

In analogy with unpolarized p-n junctionsj23 there are several regions with distinct transport characteristics in 
spin-polarized magnetic p-n junctions, as illustrated in Fig. (i) the depletion layer with space charge and large 
carrier and spin drift and diffusion; (ii) the carrier diffusion regions (CDR) which are neutral and where the minority 
carriers' drift can be neglected. CDR are characterized by carrier diffusion lengths L np for electrons on the p side and 
L pn for holes on the n side; (iii) the spin diffusion regions (SDR), which are neutral and where spin (both majority 
and minority) drift can be neglected. SDR are characterized by spin diffusion length L sp on the p side and L sn on 
the n side; (iv) the homogeneous regions in the rest of the junction, which are neutral, and where the carrier and spin 
densities assume their equilibrium values. There is no diffusion, only the majority carriers' drift. 

This section presents a unified picture of carrier and spin transport in magnetic p-n junctions. We first describe the 
profiles of carrier and spin densities inside the bulk regions, as dependent on the densities at the depletion layer, which 
are calculated next by modifying Shockley's model to the spin polarized case. The four (not independent) important 
assumptions used are (a) neutrality of the bulk regions, (b) small injection of the carriers across the depletion layer 
and at the biasing contacts, (c) the existence of a thermal quasi-equilibrium across the depletion layer even under 
applied bias and source spin, and (d) continuity of spin current across the depletion layer. Our analytical results, 
summarized in Table [n| show how the carrier density is influenced by both bias (as in the unpolarized case) and 
nonequilibrium spin, and, vice versa, how nonequilibrium spin is influenced by both bias and nonequilibrium carrier 
density. This interplay is imprinted most significantly in the dependence of the I-V characteristics of the magnetic 
diodes on nonequilibrium spin. 



A. Carrier and spin profiles 



1. p region 



In the p region the hole density is uniform, p = N a . Electrons are the minority carriers whose diffusion is governed 
by the equation 

Sn" = *L, (5) 

np 



where the electron diffusion length is L np = y 1 D np T np . We remind that if two subscripts are used in a label, the 
first denotes the carrier type (p or n) or spin (s), and the second the region or the boundary (p, n, L, or R); if only 
one subscript is used, it denotes the region or the boundary. Equation (^) is obtained by combining Eqs. (|l|) and 
(|J), neglecting the electric drift force (magnetic drift vanishes in the bulk regions), and defining l/r np = rN a . The 
boundary conditions for the electron density are Sn p at x = — w p and 5ul (yet unknown) at x = —d p . The boundary 
position of the depletion layer is not fixed, but changes with the applied voltage and the equilibrium magnetization 
(through Vb, see Appendix [A]) as 



P V q N a N a + N d W 

It is useful to introduce w p = w p — d p to describe the effective width of the p region. The solution of Eq. (||) can then 
be written as 



where i] np = (x + d p )/L np and 



Sn = 6n L cosh(f?„ p ) + F np sinh(r?„ p ), (7) 

_ Sn L cosh(w p /L np ) - Sn p 

£np — . r 7 ~ 77 — x ■ (°) 

smh(w p /L np ) 
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"Flux" parameters F are central to our analysis, since they determine the currents at the depletion layer. Effectively, 
F measures the change in the nonequilibrium (here carrier Sn) density over the length scales of the (here carrier L np ) 
diffusion length: For a short p region, L np w p , F np w (Sul — Sn p )L np /w p , while for L np -c W p , F np w Sul- The 
electron current profile, J„ = —D np Sn', is 



J n = [5n L sinh(77„ p ) + F np cosh(?7„p)] . (9) 



^ n p 

At the depletion layer, x = ~d p , the current is 

D 



JnL — —-j—Fnp- (10) 

The spin density is also described by a diffusion equation. From Eqs. (g) and (Q), under the conditions of charge 
neutrality and magnetic uniformity, we obtain 

6s =yt ~ a °P TT' ( U ) 

^sp ^lp 



where L\ p — D np T\ p and the effective spin diffusion length in the p region is L sp — D np T sp , where l/r sp = 
l/^np + i s the effective spin relaxation rate, reflecting the fact that, in addition to intrinsic spin relaxation 

processes, carrier recombination degrades spin. The second term in the RHS of Eq. ( |TT| ) acts as a local spin source, 
and appears because a change in the electron density, Sn, drives spin by intrinsic spin relaxation processes to ao p Sn 
[see Eq. (Q)], thereby preserving the equilibrium spin polarization, but not the spin itself. The boundary conditions 
for the spin density are 5s p = 5s(—w p ) and, yet unknown, 6sl = Ss(—d p ) . The solution of Eq. ( [TT| ) is 

5s = 5s L cosh(r? sp ) + F sp sinh(?/ sp ) + a 0p 5n, (12) 

where rj sp = (x + d p )/L sp , Ssl = Ssl — ao^SriL is the effective nonequilibrium spin at L, and 

Ss L cosh(w p /L sp ) - Ss p 
fsp - ■ 7 } ~ it — s {^) 

is a normalized spin flux with Ss p = Ss p —ao p Sn p . For a large spin diffusion length, L sp ^> w p , F sp m (Ssl—5s p )L sp /w p , 
while for L sp w p , F sp w 6s l- The first two terms in the RHS of Eq. (|l2| ) describe the deviation of the spin density 
from ao p n, while the last term represents the deviation ao p 6n which is solely due to intrinsic spin relaxation (Xi) 
processes. The spin current, J s = —D np 6s', has the profile 

J s = - El!f- [Ss L sinh^p) + F sp cosh(r/ sp )] + a 0p J n . (14) 

Lisp 

The first two contributions describe the spin flow due to spatial variations in 5s, while the last term represents the 
spin flow associated with the spin-polarized electron current. Finally, at the depletion layer, x = —d p , the spin current 
is 

JsL = T^-Fsp + OLQhJnL- (15) 

l->sp 

The first term can be neglected if the spin polarization is close to its equilibrium value (which is typically th e cas e at 



small biases and no source spin). The second term is important for spin extraction at large biases (see Sec. IV E ) 



2. n region 

In the n region only spin diffusion needs to be examined, as to a very good approximation n = Nd (charge neutrality 
actually requires that n = Nd + 6p, where 8p is the deviation of the hole density from equilibrium; this gives a small 



contribution to spin density 5sr, as is discussed in Sec. IV E and Appendix Ch. Electron spin diffusion is described 



by the equation (obtained from Eqs. (g) and (j^) neglecting electric and magnetic drifts and recombination processes, 
as p < N d ) 



(16) 
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where L sn = y/D nn Ti n . We introduce w n — w n — d n as the effective width of the neutral region, with bias and 
equilibrium spin polarization dependent depiction layer boundary 



2eN a V b -V 

dn = ^W d N^N- d (17) 
The boundary conditions for the spin density are Ssr = Ss(w n ) and 5s n = Ss(d n ). The solution of Eq. ( pif ) is 

5s = 5sRCOsh(r) sn ) + F sn sinh(?7 s „), (18) 

where n an = (x - d n )/L sn and 

5s n - Ss R cosh(w n /L sn ) 
smh{w n /L sn ) 

The normalized flux is F sn ss (Ss n — 5sn)L sn /w n for a short n region, L sn 3> w n , while F sn « — Ssr when L sn <C w n . 
The spin current, J s = —D nn Ss', is 

J s = [Ss R sinh(rj sn ) + F sn cosh(?7 sn )] . (20) 



Finally, at the depletion layer, x = d n , the spin current is 

J s r — — f nn F sn . (21) 

The spin current at the depletion layer boundary is solely the diffusion current due to a spatially inhomogeneous 
nonequilibrium spin in the region. Electrons with just the equilibrium spin polarization will not contribute to spin 



flow within the model approximations (see Sec. IV E for a discussion of how the neglected terms affect the carrier and 
spin transport). 

B. Carrier and spin densities at the depletion layer 

Let <fr(x) be the electrostatic potential resulting from the application of applied bias V (that is, not including the 
equilibrium built-in potential V&). We assume that all the applied bias drops within the high resistance, carriers 
devoid, depletion layer: 

0(d„) - <f>(-dp) = V, (22) 

so that (j) is constant in the bulk regions. Further, let fi be the deviation of the nonequilibrium chemical potential 
from its equilibrium value; [i is generally spin dependent: we will denote it as \i\ for spin up and /X| for spin down 
electrons. That /i is a good description of the carrier and spin off-equilibrium energy distribution follows form the 
well established fact that energy and momentum relaxation proceeds much faster than carrier recombination and spin 
relaxation. For a nondegenerate statistics, spin up and spin down electron densities can be written as 

n^(x) = nio(x)exp[<j)(x) + fj,i(x)], (23) 
n l( x ) = n l0 (x)exp[<j)(x) +m(x)], (24) 

where rifo and are the equilibrium values; we have made explicit the fact that all the quantities describing the 
densities vary in space. The electron, n = n-j- + and spin, s = n-f — nj_, densities are 

n = exp(0 + jj, + ) [no cosh(/i_) + sq sinh(/i_)] , (25) 
s = exp(</> + [no sinh(/i_) + sq cosh(/i_)] , (26) 

where fi± = (fx-\ ± n{)/2. Finally, the spin polarization 

tanh(/x_) + otQ 



1 + a$ tanh(^_) 



(27) 



depends on /i_ only (while n and s depend on both /i + and U-). 
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Substituting Eqs. (J2j^) and (|26|) into the equations (|lj) and (|2j) for the electron carrier and spin currents, we obtain, 

J n = -D n {npb + + sfi'_) , (28) 
J s = -D n (V- + W+) ■ (29) 

It may be tempting to associate /i+ with only charge, and /i_ with only spin (as done, for example, in Ref. |9l| ). It 
would then follow from Eq. ( p8| ) that in a semiconductor with a uniform carrier density a charge current would flow (or 
a spin emf would appear) if a nonequilibrium spin gradient (or, equivalently here, spin polarization gradient) would 
be maintained.tll This is wrong, as can be se en dir ectly from Eq. (uh which shows that spin can contribute to charge 



current only through magnetic drift (see Sec. IV F ), Although /i_ indeed suffices to determine a, it also influences 
n. If n is to be uniform and //_ has a finite gradient, then /i + must change to ensure that n is unchanged, as follows 
from Ea_ (p5|). However, a spin emf due to spin polarization gradient would appear in degenerate semiconductors or 
metals,E2l as mobilities and diffusivities for spin up and down species would generally be different in this case, and 
spin diffusion directly affects charge current .O 

1. Shockley's condition of constant chemical potentials 

We now apply the condition of constant chemical potentials in the depletion layer to connect the charge and spin 
densities at the left (L) and right (R) depletion layer edges. First notice that 

, , . a(x) — ao(x) , s 

tanh(u_J = — -hV = const, (30) 

1 — a{x)ao(x) 

from which follows that the spin polarizations at L and R are connected without an explicit dependence on bias. We 
will now express ol, til, and sl in terms of the nonequilibrium spin polarization in the n region, SaR = Ssr/N^; we 
will evaluate SaR explicitly from the input parameters in the next section. 
It follows from Eq. <M) that 



a 0L (l - a 0R ) + Sa R (l - a 0L a 0R ) 

aL = 1 2 — ri 7 ^ ' y 61 > 

1 - a 0R + ° a R\ a 0L - a-OR) 



If Scvr — 0, then ah = cvql- In other words, only nonequilibrium spin can be injected from the majority region through 
the depletion layer. In the case of a homogeneous spin splitting (aoL — cton)i $ a L = Sa R , that is, the nonequilibrium 
spin polarization is constant across the depletion layer. Also note that aL depends on the applied bias only implicitly, 
through the possible bias dependence of da r. 

The carrier and spin densities at L are determined by both 5a R and V. Equations ( ^5|) and (|2^) yield 

v fi , x a ° L ~ a ofl.\ / QO x 
n L = n 0L e 1 + ba R - 2 — J, (32) 



sl = s 0L e v 1 H j • ( 33 ) 

V a 0L 1- a^ R ) 

In the absence of nonequilibrium spin (San = 0), the, above formulas reduce to the well known Shockley relation for 
the minority carrier density at the depletion layer ,c3 = noiexp(V r ), and the analogous formula for spin, sl — 
soLexp(I^), so that the equilibrium spin polarization aL — aoL is preserved. Equations ( |3^ ) and (^) demonstrate 
the interplay between charge and spin in magnetic p-n junctions: nonequilibrium spin 5a R can significantly affect the 



minority carrier density (thus the junction I-V characteristics, as will be shown in Sec. Ill C ) and spin, while bias affects 
both the carrier and spin densities. If the band spin splitting is homogeneous («ol = ^or)^ nonequilibrium spin does 
not influence the minority carrier density [and affects the spin density in a trivial way: sl — noL&L exp(V)]. Equation 
( |32| ) suggests that the charge response, Sul, to nonequilibrium spin can be maximized by maximizing the difference 
in the equilibrium spin polarizations, \aoL — «o_r|, and having aoR as close to ±1 as possible (the case of aoR = ±1 
is pathological, and is excluded from our theory by the assumption of small injection, whereby ul N a , Nj). 

2. Continuity of spin current in the depletion layer 

In the previous sections 5a r was treated as an unknown input parameter to obtain the carrier and spin profiles, 
and specifically the carrier and spin densities at x = —d p . Calculation of SaR is performed in this section. The 
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knowledge of 5ctR will complete the formalism necessary to calculate any quantity of the magnetic p-n junction under 
general conditions of applied bias and source spin, with the stated constraints of the model. In the spin-equilibrium 
case (5a = 0) the calculation made in the preceding section suffices to get all the necessary boundary conditions. The 
reason is that the carrier density in the majority side is uniform, n = N#. Spin, however, does not behave similarly to 
the majority carriers even in the majority region. Spin can be injected into the majority region, and diffuses, rather 
than drifts, there. This is why the unknown San needs to be specified by another condition. Here we apply the 
condition of the continuity of spin current in the depletion layer. Physical justification for this condition is the fact 
that in the depletion layer, devoid of carriers and spin, spin relaxation, proportional to the spin density, is inhibited. 
One can write from Eq. (Q) 

JsB. = JsL — Js, relax, (34) 

where J s pH a x is the spin relaxation current (similar to the carrier recombination current used in treating unpolarized 
junctionstSl), 



Js, relax / dx 



"(sp - s po) + - 



(35) 



We neglect J s ro i ax in the following treatment .E3 

Equations @, @, (H), and ©, together with Eq. @, for m a full, self-consistent set of equations needed to 
extract Ssr (or, equivalently, San), and thus complete the structure of the model. In the process of extracting 5sr, 
we apply the condition of low injection, and neglect the terms of the order of exp(F) when compared to Nj- The 
result is 

5s R = 7o<5s„ + jiSs p + i20iQ L 8n p - j 3 s l (e v - l) , (36) 

where the geometric/transport factors are 

7o = l/cosh(u) n /L sn ), (37) 
D np \ ( L sn \ tanh(w„/L s „) 



71 - i^iif u£ ■ (38) 

Dnp\ ( L sn ^ tanh(u;„/L sn ) 



\D nn J \L np J smh(wp/L np ) 
73 = 72 cosh(w p /L np ). (40) 

Equation (|3^) expresses Ssr in terms of the known input parameters, and can be used as an input for determining the 
carrier and spin densities at the depletion layer, as well as the carrier and spin profiles in the bulk regions. The first 
contribution to Ssr comes from the source spin at the right contact, Ss n . The second and the third terms in the RHS 
of Eq. ( p6| ) come from the source spin and the carrier densities at the left contact, and a result of spin injection by 
the minority electrons through the depletion layer. Finally, the last term, which usually is negligible, results from the 
spin flow of the minority electrons having the equilibrium spin polarization (that is, as if no spin o r minority electron 



source were present). This term, for large forward biases, leads to spin extraction (see Sec. IVE). In most practical 
cases the source spin is injected either in the majority or in the minority regions, not both. Then the contributions 
to Ssr can be considered separately, with either the first, or the second and the third terms in the RHS of Eq. 
contributing. The last term (that with 73) can be usually neglected in the low injection limit. Implications of Eq. (p£ 



for spin-polarized transport in magnetic p-n junctions will be explored in Sec. IV 
The content of this and the previous sections is summarized in Table O. 



C. I-V characteristics 



Charge current in a magnetic p-n junction is driven by both external bias and source spin. Neglecting carrier 
recombination in the depletion layer, the charge electron current is the current that appears at the depletion layer in 
the minority side, x = —d p : j n — —qJ n L- Equation ( |lC)| ) gives 



]n = hn + j In + J2n, 



(41) 
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p region carrier density and current 




6n" = Sn/L^p 




6n = Sn L cosh(ry np ) + F np sinh(?7 nf) ) 








J n = -{D n p/L np )[5n L sinh(77„ p ) + F np cosh(r; n p)] 


TJnp 


(x + d p )/L np 




[SnL cosh(w p / L„ p ) — 6n p ]/ smh(w p /L np ) 


6n L 


n 0L e v [l + Sa R (a 0L - a 0R )/ (1 - Oo«)] - n 0L 


JnL 


' (_-Dnp/ Lnp^j-Fnp 


p region spin density and current 




Ss" = Ss/Llp + 65/ Lfp 




6s = 6s L cosh(?7 S p) + F sp sinh(r; S p) + a 0p 6n 




Js — Dnp^s 


Js ~- 


= ~{D np / L S p) [6§l sinh(r; S p) + F sp cosh^p)] + oto p J„ 


T) S p 


(x + d v )/Lsp 


5§ L 


6sl — aoL6nL 


ds p 


6s p — ao p 6n p 


Fsp 


[6§l cosh(w p /Lsp) — 6s p ] / smh(vjp/ L sp ) 


6s L 


SoLe v [l + (6a R /a 0L )(l ~ a LCe OR )/(l — ct? )R )] — s 0L 


JsL 


— (D np /L sp )F sp + olqlJixl 


n region spin density and current 




6s" = 5s/L 2 an 




5s = 6s R cosh(77 sn ) + F an sinh(?7 sn ) 




Js — DnnS S 




Js = -{D nn /Ls n ) [5s R sinh(j7 sn ) + F an cosh(r; sn )] 


T]sn 


{x d n ) / L S n 


TP 

± sn 


[6s n — 6s R cosh(w n /L sn )] / smh(w„/L sn ) 


6s R 


7o<5«n + li5sp + r y2ceoL5n p — 73S0L (e v — l) 


To 


1/ cosh(w„/L sn ) 


7i 


(Dnp / Dnn ) (Lsn / Lsp ) [tanh(w„ / L S n ) / Sinh(lDp / -Lsp)] 


72 


{D np / D nn )(L an /L np ) [tanh(w n / L sn ) / sinh(w p / L np )\ 


73 


72 cosh(w)p/L„p) 


JsR 


{Dnn j Lsn)F S n 



TABLE II: The carrier and spin densities and currents in the bulk regions of a magnetic p-n junction. Only electrons are spin 
polarized (spin polarization of holes is treated in Appendices [X] and g). For both the p and n regions, the diffusion equations 
and the equations for currents, as well as the explicit formulas describing the spatial profiles of the densities and currents in 
the bulk regions are given. The notation is summarized in Table pi 



where 

jon = j gn (e V ~ 1) , (42) 

Jin = J gn e OatR— j — > ( 43 ) 

1 a 0R 

— ■ 1 (AA\ 

Jin Jgn , / ~ 7T \ • V^^l 

cosh{w p /L np ) n 0L 

By jgn we denote the electron generation current (current of thermally excited electrons in the p region close to the 
depletion layerES): 

j s „ = ^n 0i coth(^y (45) 
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The generation current depends on the equilibrium magnetization through noz, (see Appendix |A|) . A magnetic p-n 
junction works as a diode when both electrodes are Ohmic (Sn p — 0), in which case j n = jo n + ji„. This current can 
be also written as 

J0n+ Jin = Jgn , (46) 

n-OL 

a notation which emphasizes the crucial role of the minority carrier density at the depletion layer for charge transport. 
Equation (S) describes the usual rectification current, which (for an Ohmic contact) is the only carrier current in 
magnetically homogeneous junctions («ol = ckor), or in junctions lacking nonequilibrium spin (SctR = 0). Once a 
nonequilibrium spin is present, and the carrier bands are inhomogeneously spin split, the current is modified by j± n , 
the spin-voltaic current, the charge current caused by nonequilibxium spin. The spi n-volt aic current does not vanish 
at zero bias, giving rise to the spin-voltaic and spin-valve effectsLa discussed in Sec. IV C| . Including the hole current 
(see Appendix [b]) , the total charge current reads 

3=3n + 3 P - (47) 

Here we consider holes to be unpolarized, so that 

Jp-J 3P (e y -l), (48) 



with 



qDp n 



Jgp = -J poRCoth(w n /L pn ) (49) 

^pn 

being the hole generation current. The hole current is affected by magnetic field only through pqr (see Appendix |a|). 
If also holes would be spin polarized, the hole current would depend on the nonequilibrium hole spin polarization, and 
would exhibit all the spin phenomena we discuss for electrons. The corresponding formulas are presented in Appendix 

For spin injection problems it is often useful to consider the spin polarization of the charge current, not only the 
density spin polarization a. The current spin polarization is defined as a,j = j s /i, where j s is the spin current 
associated with charge flow. In our case of only electrons being spin polarized, j s = —qJ s - Since j is a conserved 
quantity, the spin polarization profile is the same as the profile of the spin current, already given in the previous 
sections. As will also be demonstrated in the discussion of particular cases of interest, aj can differ significantly from 
a. Unlike for a, for example, the magnitude of aj can be greater (even much greater) than unity (if spin up and down 
electrons flow in opposite directions). The knowledge of the current spin polarization is essential particularl y in s tudies 



of spin injection, where typically one assumes that aj is conserved across the injection interface (see Sec. IV D), as a 
result of the continuity of the spin current. 

We close this section by explaining qualitatively the physics behind the spin- voltaic current ji n . Equation (|4l]) can 
be understood pather simply by considering the balance between the recombination and generation currents in the 
depletion layer .l3 In the following we put 5n p = 0, to simplify the discussion. Let C, nn and C, np denote the conduction 
band splitting in the n and p regions, respectively, as illustrated in Fig. [I]. The recombination electron current is the 
current of the majority electrons flowing from n to p. It is essentially the current of electrons with enough energy to 
cross the potential barrier in the depletion layer. This barrier is different for spin up (Vf,j = Vb + £„„ — C, n p) and spin 
down (Vbi = V, — Cnn + Cnp) electrons. Within the Boltzmann statistics the recombination current of spin up and 
down electrons, under applied bias V, is 

.M = Kn m e- v »-^+^» +v , (50) 
Jri = Kn Rl e- v ^-^ +v . (51) 

where if is a spin-independent constant. The recombination current is proportional to the number of electrons tir 
available for thermal activation over the barrier, and the thermal activation Boltzmann factor exp(— Vb + V). 

The generation currents are the electron currents (flowing from p to n) due to the minority electrons thermally 
generated in the diffusion region on the p side (Fig. [l]), and swept by the large built-in field to the n side. The 
generation currents are bias independent, and must equal the corresponding recombination currents if V = 0, so that 
no net current flows in equilibrium. Thus 

Jgl =Kn om e- v "-^+^, (52) 
j gi = Kn 0Ri e- Vb+ ^-^. (53) 
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The total electron charge current, j n = j r t — jg] + irj. — igti can ^ e expressed through the equilibrium and nonequi- 
librium electron spin polarizations, using formulas from Appendix |X| The result is 



j n = Kn 



OL 



1 



OtQL — CKQR 



1 



a 



OR 



- 1 



(54) 



which is, up to a constant, Eq. (41) (the constant K, which is proportional to the generation current, can be obtained 
rigorously only by solving the corresponding diffusion equations) . The above reasoning explains the spin- voltaic effect 
in magnetic p-n junctions as resulting from the disturbances of the balance between the generation and recombination 
currents. The nonequilibrium spin itself, 5a R , which is an input for Eq. (^), must be obtained by considering the 
full set of assumptions leading to Eq. (36). 



IV. DISCUSSION 



As an application of our theory we discuss several important manifestations of spin-polarized bipolar transport in 
magnetic p-n junctions, and illustrate the examples numerically with GaAs materials parameters. The specific cases 
we consider are spin injection (through the depletion layer) by the majority carriers, spin pumping by the minority 
carriers, the spin-voltaic effect, external (source) spin injection by the biasing electrode, spin injection and extraction 
at large biases, and magnetic drift in the neutral regions. 

The reason for choosing GaAs for numerical examples is that GaAs is the best studied semiconductor for spin 
propertiesB Spin can be injected into GaAs both optically and electrically, and high quality magnetic hybrid semi- 
conductor staiGtures based on GaAs can be potentially fabricated, as underlined by the discovery of ferromagnetic 
(Ga,Mn)As.Eil)23 For integration with semiconductor technology, however, it would be much more desirable to have 
Si-based spintronic devices. Although optical orientationel of electron spins in Si is not effective because of the band 
structure (unlike GaAs, Si is not a direct band-gap semiconductor), there seems to be no fundamental reason why-spin 
could not be injected into Si electrically; thus far, however, electrical spin injection into Si has proved elusive.Ej In 
addition to the economic reasons of easy technological integration, Si could offer other advantages over GaAs, such as 
(expected) JjOager spin relaxation times (due to the weak spin-orbit coupling and the absence of the D'yakonov-Perel' 
mechanisrnQii3 of spin relaxation in centrally symmetric Si), and much larger intrinsic carrier density rii (important 
for bipolar conduction). 

The numerical examples in the following sections are all based on a symmetric GaAs magnetic diode with the fixed 
parameters N a ~ Nd = 10 16 cm~ 3 , n, = 1.8 x 10 6 cm -3 , D n = 100 cm 2 /s, and r„ = T\ = 1 ns (equal in both 
regions), w p = w n = 3 //m. The derived parameters are L np 3.2 /im, L sp 2.2 /im, and L sn rs 3.2 /im. Other 
parameters (bias, equilibrium and nonequilibrium spin) will be specified according to the physical situation. The 
materials parameters are for room temperature, so the chemical potentials will be given in the units of ksT w 25 
meV. 



A. Spin injection by the majority carriers 

Under the low injection conditions nonequilibrium spin cannot build up in magnetic p-n junctions, as was shown in 



Sees. Ill B 1 and IIIB2. Only if a nonequilibrium (source) spin is externally injected into cither region of the junction, 
spin injection through the depletion layer is possible. Here we consider the case with a magnetic n side ((Xor ^ 0) and 
a nonmagnetic p side (ccol = 0), and inject the source spin at the right contact (but not by the contact itself), so that 
5a n 0. The left contact remains Ohmic with equilibrium carriers and spin (Sn p — Ss p = 0). The nonequilibrium 
spin at the depletion layer in the n region is obtained from Eq. (pBq ) (see also Table ||) as 

Ss n 

5s R = Sa R N d = " — r. (55) 

cosh(w„/L s „) 

This boundary condition for spin at the depletion layer cau be physically formulated by requiring that the spin current 
of the majority carriers vanishes at the depletion layer.Efl This is quite natural to assume, since the spin current in 
the n side is proportional to <x n N^ while the spin current in the p side is proportional to the much smaller a^n^. 
Since J Sj r = J s l, we can neglect J s r relative to J s in the rest of the n region. Eq. (p5| ) then follows. 

On the left side of the depletion layer the Shockley condition, according to Eqs. (|32|) and (p3[), gives for the electron 
and spin densities 

n L = n 0L e 1 - ^— , (56) 
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FIG. 2: An example of the majority carrier spin injection through the depletion layer (shaded). Shown are the spatial profiles 
of the electron (solid) and spin (dashed) densities in the magnetic p-n junction described in the text, with aofl = 0.5, 8a n = 0.4 
(the p region is nonmagnetic and the left electrode remains Ohmic), and forward bias V=+0.8 volt (~ 32 kgT). The left vertical 
axis is for the densities, while the right axis is for the spin polarization, which is represented by the dotted lines labeled with a. 
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(57) 



(58) 



If the source spin has the same direction of polarization as the equilibrium spin in the n region, the electron density 
ri£, and thus the current through the junction, is reduced. If they are antiparallel, and the current are enhanced. 
Neither spin polarization, San nor cvl, depends on V (except for a small dependence through w n ), being the same 
for forward and reverse biases. In nonmagnetic junctions (aoa = 0), all the nonequilibrium spin polarization is 
transferred to the minority region, aj, = <5cv_r, where the nonequilibrium spin has no effect on charge and current, 
since Sul = V&l exp(V). This case has been studied numerically for a realistic model of a spin-polarized nonmagnetic 
p-n junction.o 

The reason for the absence of spin injection through the depletion layer from a magnetic n region to the nonmagnetic 
p region, without a source spin, is the balance between the carrier densities and thermally activated processes of forward 
conduction. Let the n region be positively magnetized, so that there are more spin up than spin down electrons. For 
a nondegenerate statistics, the number of spin up (down) electrons depends on the spin splitting (2£ nn ) of the band 
as exp(£„„) [exp(— („„)]. In the forward transport, electrons need to be thermally activated to cross the barrier of the 
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FIG. 3: Calculated current spin polarization for the majority carriers spin injection. The same parameters as in Fig. g 
apply. Both the current spin polarization aj and the density spin polarization a are shown for comparison. The current spin 
polarization is enormous in the n region, decreasing upon reaching the depletion layer, and staying smaller than a in the p 
region. 



built-in voltage lowered by the external bias. The barrier height is, however, different for spin up and down electrons. 
Indeed, spin up (down) electrons have the barrier higher (lower) by leading to the modulation of the transport 
rate by exp(— £ rm ) [exp(£ nn )]. These exponential factors exactly balance the modulation of the carrier densities. As a 
result, there is no difference between the transfer rates (density times the thermal activation probability) for the spin 
up and spin down carriers, the spin up and spin down currents are equal, and there is no spin current at R (and, by 
the continuity of the spin current also at L) and thus no spin injection into the minority region. 

Figure ^| shows the electron and spin densities, using our model equations (Table |n|) , for the GaAs magnetic junction 
example, with aoj? = 0.5 and Sa n — 0.4, and a forward bias of +0.8 volts. Spin injection into the minority region 
is very effective; qj, is slightly greater than dan [due to the denominator in Eq. (|3lJ)]. A comparison between the 
current spin polarization (the profile is the same as for the spin current J s ) and the density spin polarization is in 
Fig. [| The current polarization is huge at the point of spin injection, since in order to reproduce the spin polarization 
a n by electrical spin injection (which would depend essentially on aj, see Sec. IV D), aj would need to be that large. 
This is of course not possible, since electrical spin injection from a ferromagnetic electrode provides aj < 1, since aj 
in ferromagnets is close to the density polarization there. The current polarization decreases upon approaching the 
depletion layer, since there the spin current decreases in order to be equal to the spin current at L, which is driven by 
the much smaller density of the minority electrons. Figure |^ shows the chemical potential profiles for the case. The 
chemical potentials are chosen to be zero (similarly to <ft) at x = w n in the spin unpolarized (but biased) junction, so 
that at x = w p they are — V if the contact is Ohmic, as is the present case. This is the cause of the rapid decrease 
of the /i's to fi — —V at x = —w p . Spin injection in this graph is visible from the finite value of /i_ (which becomes 
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FIG. 4: Calculated chemicalpotential profiles in a magnetic p-n junction under the majority carrier spin injection regime. The 
same parameters as in Fig. pj apply. The chemical potentials are expressed in the units of ksT. 

zero only in the very proximity of the left contact) in the p region. 



ic field to the majority side, 
~ ,3 ' since the spin accumulation 



B. Spin pumping by the minority carriers 

If large (source) spin density is externally injected along with the carrier density into the minority region, the 
nonequilibrium spin can reach the depletion layer and be swept by the built-in e 
where it accumulates. We have named this effect minority electron spin pumping, 
(which is also a spin amplification, considering that the resulting spin in the majority region is much larger than that 
in the minority region) depends on the intensity of the spin current of the minority carriers. The faster the carriers 
arrive at the depletion layer, the mare spin accumulates in the n side. In effect, this is an analogue of the optical 
spin pumping in the majority regionJJ except that the role of circularly polarized light is played by the spin polarized 
minority carriers. 

As an illustration consider a nonmagnetic spin-polarized p-n junction (aoz, = ao-R = 0). Let the carrier and spin 
densities at the left electrode only be out of equilibrium: 5n p , Ss p ^ 0. This -happens, for example, when a junction 
is illuminated by circularly-polarized light (like in a spin-polarized solar celEa) or if the junction is part of a spin- 
polarized junction transistor, in which case the left electrode simulates the action of the emitter. Equation (|36| ) gives 




FIG. 5: An example of a minority carrier spin pumping through the depletion layer. The junction is nonmagnetic, but spin- 
polarized, and the carrier and spin source is placed at the left electrode, giving 5n p = 1 x 10 15 cm' 3 and Ss p = 1 x 0.9 x 10 15 
cm -3 (Sa p = 0.9). A reverse bias of —0.8 volts is applied (increasing the width of the depletion layer compared to Fig. ^). 



the "pumped" spin polarization in the majority side as 



D n p\ ( Lsn\ tanh(w„/i s „) 



\DnnJ \L S p J smh(w p /L sp ) 

For a large majority region, w n ^> L sn , the injected spin is (below only holds if w p << L sp ) 6s r m 
[Pnp/ D nn )(L sn /wp)5sp, while for a short majority region, w n -C L sn , the injected spin is 6s R m 
{D n p/D nn )(w n /Wp)5sp. The amount of the pumped spin polarization, relative to the amount of the source po- 
larization is 

6a R ( D np \ f L sn \ ta,nh(w n /L s 



(61) 



,D nn J \L sp J smh(wp/L sp ) N d 
Spin pumping is most effective when the p region is short, w p -C L sp , when 

6a R _ ( D np \ mm(L sn ,w n ) 6n 



Dnn J Wp N d 



If both L sn and w n are significantly greater than w pi the pumped spin (and even the spin polarization) can be 
comparable to the source spin (source spin polarization). 
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FIG. 6: Calculated current spin polarization in the minority spin pumping regime. Both current, aj, and density, a, spin 
polarization profiles are shown. The current spin polarization starts at a value larger than 1 at x = —w p , remains constant 
across the depletion layer where the spin current continuity is assumed, and decays somewhat in the n region, where its 
magnitude is much larger than that of density spin polarization. 



A qualitative argument for the spin pumping is as follows. In the minority (p) side, the spin current goes roughly 
as D np 5sp/wp, where we chose the largest spin in the region (being the source spin Ss p ) and the smallest length scale 
for the spin decay (here w p ). On the n side the spin current, along similar reasoning, would be s» D nn 6sn/ L sn , where 
Ssr is the largest spin in the region and we chose L sn to be the smallest length scale. Equating the two currents gives 
Eq. (|6l|). Put in words, spin carried by the minority carriers arriving at the depletion layer is swept into the majority 
region by the large built-in field. In the majority region the spin both diffuses away and relaxes. In a steady state, 
the incoming spin flux must equal the outcoming diffusion and relaxation, which are proportional to the spin density, 
so that the greater is the spin influx, the greater the spin density. 

A numerical example is shown in Fig. |^. The source carrier and spin densities are 5n p — I0 15 cm" 3 and Ss p — 
0.9 x I0 15 cm" 3 (the spin polarization a p — 0.9). The junction is under reverse bias of —0.8 volts (note the increase 
width of the depletion layer compared to Fig. |2|). The pumped spin polarization ctR is about 5%. In our numerical 
example all the length scales involved are comparable (roughly 3 /an), diffusivities uniform (D np — D nn ) so Ssr rs Ss p . 
In Fig. ^ we plot the current spin polarization aj to demonstrate that it significantly differs from the density spin 
polarization a. In this example aj is larger than I at the left electrode due to the chosen boundary conditions, 
and in the n region it is much greater than the density spin polarization. The chemical potential profiles for the 
case are shown in Fig. [?]. In the majority region /i + nearly vanishes, while /if rs /tj, demonstrating a positive net 
nonequilibrium spin polarization in the n region. The small magnitudes of the noncquilibrium chemical potentials in 
the majority region still yield large spin density, since they appear in the exponent which multiplies the equilibrium 
carrier density, which is large in the majority region (and small in the minority, where the chemical potentials have 
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FIG. 7: Calculated chemical potential profiles in a nonmagnetic spin-polarized p-n junction under the minority carrier spin 
pumping regime. The parameters as in Fig. |E| apply. The input shows the majority region values on a scale where different fi's 
are visible. 

accordingly large magnitudes). 



C. The spin- voltaic effect 



A spin- voltaic effect is a generation of charge emf or current 
voltaic effect was the Silsbee- Johnson spin-charge couplingl 



jy noncquilibrium spin. A first realization of the spin- 
in a ferromagnetic/nonmagnetic metal interface with 



nonequilibrium spin injected into the nonmagnetic metal. The emf across the interface arises due to the difference in 
the chemical potentials in the two metals, with different effects on the different spin states. Analogous phenomena 
can occur in many other hybrid systems (semiconductor/metal or semiconductor/semiconductor). Here we describe 
a specific realization of the spin-voltaic effect in magnetic p-n junctions, where the role of the interface is played by 
the depletion layer. 

Consider a magnetic/nonmagnetic p-n junction, with the p region magnetic (Sn p = 0) and the n region nonmagetic 
but spin polarized (Sa n ^ 0). No external bias is applied (V = 0). ft follows from Eq. (^) that 



da i 



Sa r 



cosh(w n /L sn ) ' 



(62) 
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which is the same as Eq. ( |55| ) (simply expressing the fact that the polarization is bias independent). As a result, there 
will be nonequilibrium carrier and spin densities in the minority region [see Eq. ( |32"|) and ( 133 " 



5n L = n 0L a LSaR, (63) 
Ss L = n 0L 8a R . (64) 

The nonequilibrium minority carrier density Shl leads to the minority diffusion and relaxation, and thus to the charge 
current (or voltage in an open circuit). The spin-voltaic current is [see Eq. (|43"|) 1 

in = j gn aoLSa R . (65) 

The current is of the order of the generation current, and changes sign with reversing either the magnetic field (which 
reverses «ol) or the orientation of the source spin Sa n . Neglecting the variation of 5a R with bias (through w n ), the 
open circuit voltage for the spin-voltaic effect is obtained by requiring that j vanishes: 

-a 0L Sa R ) . (66) 



Jgn + Jgp 



The voltage, which is typically of the order of ksT, is negative (reverse biasing) if the polarizations are parallel, 
and positive (forward biasing) if they are antiparallel. The spin-voltaic effect here is similar to the photovoltaic 
effect, where the photocarriers generated within the carrier diffusion length L np of the depletion layer are swept by 
the built-in field to the majority side, generating photocurrent. A spin-voltaic effect arises if nonequilibrium spin is 
generated within the spin diffusion length L sn of the depletion layer, disturbing the balance between the generation 
and recombination currents. 

Indeed, in equilibrium both the generation and the recombination currents in a magnetic p-n junction are equal 
and there is no net charge flow. Let aoi be positive. Then the barrier for the majority electrons to cross the depletion 
layer (see Fig. |l|) is smaller for spin up than for spin down electrons. If the spin in the majority region is driven off 
balance (without applying an external bias), than the delicate balance of the generation and recombination currents 
is disturbed, resulting in a net charge current. Increasing the number of spin up majority electrons, for example, 
increases the recombination current, since more electrons have now a smaller barrier to cross (the generation current 
does not depend on Sa R or bias). In our geometry, the net electron flow is forward (from the right to the left, j n > 0). 
If, on the other hand, we increase the number of spin down electrons, more electrons have now a higher barrier to 
cross, reducing the recombination current, resulting in a net reverse flow (from the left to-the right, j n < 0). The 
spin-voltaic effect is the reason for the giant magnetoresistance of magnetic p-n junctions,Lj since when a bias V is 
applied, the spin-voltaic current grows as exp(V_), similarly to the normal rectification current. 

The spin-voltaic effect is illustrated in Figs. |q-[iT[ First consider parallel spin polarizations, aoL = Sa n — +0.9. 
There is no bias, V = 0. The carrier and spin densities and the spin polarization are plotted in Fig. ^. The induced 
nonequilibrium spin and charge in the p region are greater than the equilibrium values, leading to a forward current 
of electrons. The spin polarization is also higher than in equilibrium. The chemical potential profiles are shown in 
Fig. |^. If the spin polarization of the source spin is reversed, aoL = —Sa n = +0.9, the carrier and spin densities and 
the spin polarization decrease in the minority region, leading to a reverse electron current. The density profiles for 
this case are in Fig. [h], and the chemical potentials are plotted in Fig. [n|. 

Let the magnetic field B controls the conduction band spin splitting. Then ao(B) = —ao(—B). Keeping 5a R as 
an independent (of B) parameter, the direction reversal of the magnetic field results in a change in charge current: 

UB)-U-B)=j J nL{B) - 5nL{ - B) . (67) 

This is a realization of giant magnetoresistance (GMR) in magnetic diodes. The relative change of the charge 
current upon reversing the direction of magnetic field (while keeping 6a R unchanged) can be measured by the giant 
magnetoresistance parameter, here denoted as f3: 

,. ^ ( ff-*> , 

which at forward bias and exp(V^) 3> 1, in terms of the known parameters, can be expressed as 

(3 = 2Sa R - g — — -. (69) 

1 - a 0R + ° a R{ a 0L ~ a 0R ) 




FIG. 8: The spin-voltaic effect in a spin-polarized magnetic p-n junction. Shown is a junction with a magnetic p region 
{ctoL 7^ 0) and a nonmagnetic n region (oor = 0). No bias is applied. Both electrodes are Ohmic, except that there is a spin 
source at x — w„. In the example ool = +0.9 and Sa n = +0.9. The carrier and spin densities in the p region are very close to 
the equilibrium values, with a small variation due to the nonequilibrium spin. The inset shows this variation on a 10000 times 
increased scale. Both densities are higher than in equilibrium, leading to a forward charge current. 



The GMR effect is possibfe only in magnetically inhomogeneous p-n junctions with nonequilibrium spin. As a special 
case consider the p region magnetic (aofl — 0)- Then 



OL 



1 + 8a R a L 



(70) 



This case is a semiconductor analogue of the Silsbee- Johnson spin-charge coupling,E3 where a spin emf arises from 
the proximity of a nonequilibrium spin in a metal and a ferromagnetic electrode. A detailed numerical study of the 
GMR effect in magnetic diodes can be found in Ref. ^4]. Putting reasonable parameters <5«ol = 0.5 = — Scxr, Eq. J70| ) 
gives j3 = 2/3, which is a 66 % GMR. A more optimistic set, ool = 0.9 = — Solr, leads to f3 ~ 8.5, or a 850 % GMR, 
demonstrating a great technological potential of magnetic p-n diodes. 



D. Spin injection by the biasing electrode 

Thus far s n was a free input parameter of the model. If, however, the biasing electrodes themselves can inject 
spin (for example if they are magnetic) , then the source spin density will not be a good starting boundary condition. 
We consider an example of the source spin injection by the right electrode into the nonmagnetic majority, n, region, 
keeping only the p region magnetic. We assume the model in which the spin current across the electrode/n- region 
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FIG. 9: Calculated chemical potential profiles in a spin-polarized magnetic p-n junction under the conditions specified in Fig. ^. 



interface is preserved. In this scenario the boundary condition at x — w n reads (all the current at the contact is 
carried by electrons, since the hole density is in equilibrium there) 

(71) 



3 s 



where j sn = —qj s (w n ) and aj n = aj(w n ) is the spin injection efficiency (here the current spin polarization at the 
contact) equal, in an ideal case, to the spin polarization in the electrode material reduced by interface spin relaxation. 
Our strategy is to convert this boundary condition to the condition on the spin density: We calculate Ss n as a function 
of aj n j and then use the formulas derived earlier to obtain the charge current in a self-consistent manner (this is 
needed because the boundary spin depends on the current which, in turn, is calculated using the boundary spin). 
Equation (EG) gives 



3s 



qD nn \ Ss n cosh(w n /L sn ) - Ss R 



(72) 



smh(w n /L sn ) 

If we further assume that the left contact is ohmic, by substituting Eq. ( |55| ) for 5sr, the above equation can be solved 
for the source spin density with the result 

(73) 



OtJnjLsn 

ds n = coth 

qD nn \ L, 



72 sol (e 



tanh 2 (w n /L sn ) 
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FIG. 10: The spin-voltaic effect in a spin-polarized magnetic p-n junction. The same conditions as in Fig. g apply, but the 
direction of the source spin is reversed, Sa n = —0.9. The figure shows the negative spin density (— s) in the n region (and 
normal in the p side). The carrier and spin densities have values close to the equilibrium ones, but are now somewhat smaller, 
due to the presence of the antiparallel nonequilibrium spin. This density variation, which is seen in the inset on a 10000 times 
greater scale, leads to a reverse charge flow. 



The nonequilibrium spin polarization at the depletion layer then is 

„ . L sn 1 

osr = a Jn j- 



qD nn smh(w n /L sn ) 
where another geometric/transport parameter is introduced: 

73 



74 sol (e 



v 



74 



(74) 



(75) 



tanh 2 (w n /L sn ) 

In a first approximation one can put j w j 0n + j p for the current in Eq. [74]. The injected spin 5sr is then of the order 
of the minority electron density times L sn /L np . This is generally larger than the spin extraction factor coming from 
the term with 74, but still small to lead to a significant modulation of nonequilibrium spin. The charge current is 
obtained by solving Eq. p^ ) for j, with 5a r from Eq. |74|. The result is 



i ~ (if- + jon) (1 + a.7«aoi) - HOiQ L j gn (e v - l) 



N d 



(76) 



where 




FIG. 11: Calculated chemical potential profiles in a spin-polarized magnetic p-n junction. The same parameters as in Fig. p"o| 
apply. 



Unlike in the case of independent external spin source, the spin injected by the biasing electrode is very small, 
because of the small current flowing in the junction (the current is carried by the minority carriers), so that only 
a small spin current can build up the source spin. As a result, the spin injected from the contact has a very small 
effect on the charge properties of the junction. Charge current, in particular, is only slightly modified from the 
spin-equilibrium value of j = j p + j 0n . The spin- voltaic effect is absent (except for the small effect caused by spin 
extraction), since at zero bias no nonequilibrium spin is injected. Nevertheless, even if small, the TV characteristics 
modification should be observable at reasonably large biases, and could be used to characterize spin properties of the 
junction. Furthermore, the above model of spin injection, based on Eq. (|7l|), is rather simple and we cannot exclude 
the possibility of a different behavior (especially more pronounced spin dependent effects) with realistic interfaces. In 
fact, our method shows a way how to characterize spin properties of real (electrode) interfaces by measuring charge 
response of the junction. 
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E. Spin injection and extraction at large biases 



We have shown numerically in Ref. [54| that spin can be injected and extracted through the depletion layer at large 
biases, even with no source spin present. Significant spin injection from the magnetic n region into the nonmagnetic p 
region occurs at large biases and, similarly, significant spin extraction occurs from the nonmagnetic n region into the 
magnetic p region. These intrinsic spin injection phenomena have their origin in the low bias physics. Indeed, there is 
a (normally negligible) spin injection in the absence of source spins even in our theory. We have already demonstrated 
spin extraction in Eq. (|36|). If no source spin is present, then 

Ss R = -73S0L exp(V). (78) 

The nonequilibrium spin in the n region will have the sign opposite to that of the equilibrium spin in the p re- 
gion. The spin is more extracted the larger the bias is. Normally, at small biases, the injected spin polarization 
Sor — — 73S0L exp(V)/Nd is small (in the postulated low injection regime), but it shows the trend of spin extraction 
exponentially increasing with V towards the large bias regime. The reason for Eq. ( |78] ) is the continuity of the spin 
current across the depletion layer. Without any spin source, the spin current J s l will be given by the flow of electrons 
with the equilibrium spin polarization a aL [see Eq. |H])]. The same spin current must appear in the n region at R. 
For positive sql, the spin current in p is negative. In the n region, for the spin current to be also negative there must 
be a positive gradient of 5s and, since Ss n = 0, the spin at R must be negative: Ssr < 0. 

On the other hand, our theory thus far does not predict any spin injection from the magnetic n to the nonmagnetic 
p region. Indeed, in the absence of source spin, and with aoL = 0, Eq. ( fj6| ) gives Ssr = 0i— To explain the intrinsic 
spin injection observed at large biases (but also at small biases, although on a smaller scaleEa), we have to introduce 
electric field E into the picture. In fact, once the nonequilibrium spin becomes small compared to the equilibrium 
one in the majority region, even at small biases electric field cannot be neglected. We will quantify this condition 
below. There are two factors which need to be considered when introducing charge effects in spin transport in the 
bulk regions. First, we will include the electric drift force into the spin current and, second, we will explicitly account 
for charge neutrality by postulating that n = Nd + Sp (instead of what was used thus far, n = Nd). These two factors 
can be normally neglected at small biases, but here we use them to demonstrate the trends, namely the spin injection, 
which will become important at large V. 

Including the E-Held and the charge neutrality, the spin diffusion equation from Eqs. (Q) and (||) becomes 

Ss" + ESs'=-^--a 0n -^-. (79) 

sn sn 

We have neglected the nonlinear terms ~ sp, justifiably if {Ti/T pn )(5p/Nd) <C 1, which is quite generally the case at 
low injection. The above equation needs to be supplemented with the diffusion equation for holes, 

Sp" = (80) 

^pn 

In Eq. ( |79| ) the term with the first derivative comes from the electric drift, while the term proportional to Sp appears 
because of the neutrality condition n — Nd + Sp. The latter term acts as an intrinsic spin source, similarly to the 
term of Sn in the spin diffusion equation ( pT| ) for the minority electrons. The neutrality condition alsp guarantees 
that the electric field is uniform (E' = 0). Equation (|7^) has already been considered and solvecOEErEj without the 
intrinsic source term, which becomes important in bipolar transport at large biases. For completeness, we present the 
full solution to Eq. (179), as well as the spin current profile, in Appendix 0. 

The full analysis in Appendix |C| shows that at least for L pn <C L sn , the contribution from the charge neutrality 
(Sn = Sp), that is, from the hole density effects in spin transport, can be ne glected. In the opposite case, the 



contribution would lower Ssr, as can be seen easily by equating J s r in Eq. (08) to zero. The elec tric field, on the 
other hand, increases Ssr, ultimately leading to spin injection at large biases. Indeed, from Eq. (C8) one obtains for 
spin injection from the majority magnetic region to the minority nonmagnetic region in the absence of source spin, 
but at a finite bias, 

5s R » s r (L sE E) tanh [j^j ■ (81) 

To obtain E we can use the carrier current continuity across the depletion layer: J„£ = J n R, where J„l is given by 
Eq. (|l(J) and J n R, w D nn (—NdE — 5p' R ), with Sp calculated in Appendix |(]. We get 

E fa coth(w n /L pn ) + ^ UL coih{w p /L np ). (82) 
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The electric field is positive at forward biases, making 5sr and thus Ssl is of the same polarity as Sor. This explains 
the large bias spin injection (and the increase of spin polarization in the magnetic n region) observed in Ref. ^4. The 
electric field spin injection will also happen at small biases, but, because the field is very small (SriL/Nd, 8pn/Nd <C f ), 
the spin injection is negligible. However, the .E-field must be considered in cases where spin injection due to source 
spin leads to Ssr as small as Sn^. 

Equation ( j32|) yields the criterion for neglecting electric drift in spin transport in the n region (which is in contrast 
to the majority carriers, for which electric drift cannot be neglected). Indeed, one needs to compare the typical 
magnitudes of the spin drift (~ srE) with the spin diffusion {5s/L sn ), to obtain 

E{s r + Ssr) < Ss R /L sn . (83) 

For a nonmagnetic n region (s _r = 0) this is always the case, since EL sn <C 1 because of low injection (and reasonably 
assuming that the spin diffusion length is not much greater than the carrier diffusion lengths). For the magnetic region, 
the above condition is satisfied if Ssr 3> {EL sti )sqr, which roughly means that the nonequilibrium spin in the n region 
(appearing through the spin source, for example) should be greater than the nonequlibrium carrier density times the 
ratio of the spin diffusion length the the carrier diffusion length. This is well satisfied in the low injection regime, 
where the nonequlibrium carrier densities are small enough (even if L sn would be one to three orders of magnitude 
greater than the carrier diffusion lengths). However, the condition ( |S3| ) places the lower limit on the source spin to 
lead to pure spin diffusion in the n region. 

Finally, the neglected spin relaxation current J Sjrc iax also contributes to spin injection, more with increasing bias, 
since then the spin density in the depletion layer increases and with it spin relaxation. The difficulty in introducing 
Js, relax is that it depends on both bias and 5sr complicating the self-consistent process of obtaining Ssr in terms 
or the input parameters and bias. One may expect, though, that spin relaxation processes in the depletion layer 
will decrease the spin injection efficiency (that is, reduce Ssl) while allowing for larger Ssr. to balance the spin 
current in the minority region. Our numerical calculations, which take into account the effects of J s , relax, find that 
its contribution is indeed small at low biases .E3 



F. Magnetic drift in the neutral regions 



Our model and its conclusions thus far were based on magnetic p-n junctions with homogeneous magnetic doping 
in the neutral regions. The doping, and thus the band spin splitting and the equilibrium spin polarization, changed 
spatially only in the transition region. As a result, the magnetic drift force dropped from the calculations and the 
inhomogeneous magnetic doping affected the results only through the equilibrium spin densities. Here we take the 
next step and ask how would the physics of magnetic p-n junctions be affected if, additionally, the neutral regions 
were inhomogeneously doped with magnetic impurities (or, to the same effect, were homogeneously doped but placed 
in an inhomogeneous magnetic field). We will show that magnetic drift modifies both the spin injection through the 
depletion layer, and the I-V characteristics of magnetic p-n junctions. The effects of £' are qualitatively different in the 
majority and the minority regions, so we will discuss the two regions separately. Most of our discussion below applies 
equally to homogeneous (in relation to nonmagnetic doping) semiconductors with spin split majority and minority 
bands. 

Consider the majority, n, region first. In the presence of an inhomogeneous spin splitting of the conduction band, 
the electron current in the region is 

Jn = D nn {-nE + <„ - n'). (84) 

The current must vanish in equilibrium where n w Nd and s ~ a^Nd- This is only possible if a local electric field, 

E « aoC, (85) 

develops. The resulting electric drift needs to counter the magnetic drift. The existence of Eq is also warranted by 
the vanishing spin (J s ) and hole (J p ) currents. In the latter the electric field needs to balance the equilibrium hole 
diffusion ~ p' [p a is now spatially dependent through £ nn , see Eq. The field E , similarly to the built-in field 

in the depletion layer, is an equilibrium field, not an emf, as it does not lead to a net current. 

The origin of the equilibrium electric field in otherwise almost homogeneous charge situation (the majority carrier 
density is almost constant) can be qualitatively explained as follows. Take an n-type semiconductor doped inhomoge- 
neously with magnetic impurities in zero magnetic field. In equilibrium the chemical potential is constant. Switch on 



a magnetic field. At first, the chemical potential will vary with x through £ nn according to Eq. (A7). The sample will 



come to equilibrium by rearranging its charge as the electrons will move in the direction of decreasing [iQ n , resulting 
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in a constant chemical potential, but also in a space charge (see below) and a space electric potential opposing further 
electron motion. Then 

0o = -lncosh(C„„). (86) 



The equilibrium electric field is Eq — —4> , reproducing Eq. ( p5[ ) obtained from transport considerations. Electric 
potential 4>q bends both the conduction and the valence band. As for the conduction band, 4>q tends to straighten 
the lower spin band (say, the spin up band if £ nTl is positive) and steepen the upper spin band. At large magnetic 
fields the band bending of the lower spin band entirely eliminates the spatial variations of the band due to £„„ , while 
these variations are doubled in the upper spin band. The valence band too is affected. Originally constant, the band 
acquires spatial variation —q<fio to balance the equilibrium hole diffusion. 

In turn, the inhomogeneous Eq induces space charge pq, according to Poisson's equation: p = E'^eksT / q). In 
principle, both Eq and po need to be obtained self-consistently by solving for the equilibrium semiconductor densities 
taking into account Poisson's equation (this was done numerically in Ref. |54| for the transition region, where the Eq- 
like field is present due to the inhomogeneous magnetic doping). However, the induced local charge density is small 
enough to be neglected for most practical purposes (unlike the induced charge density in the depletion region) . Indeed, 
the induced changes in the carrier density come to p a /q w Nd((' nn ^D) 2 / cosh 2 ((„„), where Ad = y^efc^T '/N d q 2 ) is 
the Debye screening length in the majority region. For GaAs with e = I3e and at room temperature, the density is 
Pa/q < 2 x I0 5 (C4n[cm -1 ]) 2 cm -3 . If the magnetic splitting changes by k B T over a micron (so that (' nn ~ 10 4 /cm), 
we get po « 2 x 10 13 cm' 3 . This shows that for carrier densities greater than, say, 10 15 cm -3 the induced densities in 
the carrier concentrations can be neglected, and Eq. ( |S5| ) is a reliable estimate of Eq. In general, the space charge can 
be neglected if the band splitting varies by fc^T over the length-scales greater than Ad. This is in complete analogy 
with space charge considerations due to the usual carrier doping £3 Once Cnn^D 1, which is normally the case when 
a magnetic and a nonmagnetic semiconductor form a contact for spin injection, the space charge and its diffusion 
(~ p'o/q) cannot be neglected. Indeed, for £ changing over a 0.1 pm, the induced charge density is q x I0 15 cm~ 3 . 
Selective doping of semiconductors with magnetic impurities on spatial scales both smaller and larger than Ad can 
prove a useful tool for band structure and space charge engineering in designing new spintronic devices. 

Expanding about the equilibrium values for the densities and the electric field, the electron and spin currents in 
the n region become 

J n « -D nn (N d SE-N d 5aC nn +5p'), (87) 
J s w -D nn N d {a SE + a 5aC nn + 5a') , 



where SE describes only the electric field induced by the applied bias. Nonequilibrium charge neutrality, Sn = Sp, 
is assumed. In a homogeneous n-type semiconductor with finite £' m , maintaining a nonequilibrium spin polarization 
would lead to a spin emf according to Eq. ([S7|). For a constant 5a, for example, the spin emf is — foA£„„, where 



A( nn is the drop of Q nn across the sample. Spin injection is modified by the presence of £' nn in Eq. (88). Considering 



here only a special case of a constant (' nn and large spin polarization, |ao| ~ 1, the spin current at x = d n is modified 



from Eq. (21) to 



J s r = - [ -5s R a 0R C,' nn L sn Q + F a J , (89) 

where the new length scale for spin drift-diffusion is l/L snC = ^/{C^j4+l/L 2 sn ) (the spin polarization decay is then 
governed by two length scales, L sn ±, given by 1/L sn ± — 1/ L sn Q ± Q 1 and spin flux 

F a = fa w eX , P( " 5""f w/ f - 5a R cothK/L snC ). (90) 



To obtain 5s R , one can still use Eq. (|36j), but with L sn changed to L sn Q and coth(iu n /L s „) changed to coth(tt' n /i ;jn( ;) — 
a 0R(C' nn Lsn() /2. Similarly for the spin injection at large biases. Note that in the presence of magnetic impurities L sn 
will be greatly reduced, so that L sn ^ ps 2/£'. 

Although the inhomogeneous magnetic doping affects directly only the majority electrons, it modifies, through Eq, 
transport of the minority holes as well, and thus the I-V characteristics of the junction. The hole current becomes 

J P ~ D pn (a R(' nn 5p - 5p') , (91) 



where the first term describes drift of the nonequilibrium hole density by Eq. Together with the continuity equation 
for hole current describing electron- hole recombination, the above equation, again in the limit of a constant |ao| ~ I 
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leads to the hole current at x 



JpR 



coth 



J pnC, 



(92) 



where we introduced an effective magnetic drift length l/L pn Q = -\/(C' 2 /4 + l/L pn ) (the two length scales for the 
hole density decay are L pn ± given by 1/ L pn ± = \/L pn Q ± Cf /2). Since j p = qJ P R, magnetic drift directly affects the 
I-V characteristics of the junction by modifying the hole minority current. It is in the combination with external 
bias [6pr ~ exp(V)] that the magnetic drift generates current. This effect could be used in electronic detection of 
magnetic field gradients. 

Now we turn to the minority, p, region. Since the minority electron density can easily accommodate to spatial 
changes in Q np , no equilibrium electric field is needed to balance the magnetic drift force. The carrier and spin 
currents vanish at E = for the equilibrium electron and spin densities, unlike in the n region considered above. 
From Eqs. (]l])-(0), the drift-diffusion equations for the minority electrons and spin in the p region are obtained as 



^ + En< - C np s' ~ C P * = 
s" + Es' - ('n' - (In = 



T2 ■■ 

np 

5s 

T2 

np 



T 2 
^lp 



(93) 
(94) 



Transport of minority carriers is thus coupled with the transport of spin. As a result, the electron current (and 
thus the I-V characteristics) will depend explicitly on nonequilibrium spin and, similarly, spin current will depend 
explicitly on nonequilibrium charge. Below we solve Eqs. ( [93] ) and (^4|) for the specific model of a linear £ np (that 
is, (' np = constant) and in two limits of slow and fast spin relaxation. We will also neglect the electric field which 
is by about 5riL/Nd <C 1 smaller than the inverse of the typical decay length of the densities. Magnetic drift brings 
a new length scale, L np ^, given by \jL np ^ — \/XCnp/^ + V-^np)- The density profiles then decay with two length 
scales, L np ±, which arc the inverse of \/ L np Q ± C« P /2, depending on whether the diffusion is parallel (minus sign) or 
antiparallel (plus sign) to magnetic drift. 

We now consider the limit of vanishing 1/Li p , which corresponds to slow spin relaxation (spin diffusion length is 
the largest length scale in the problem). We will not present the full density profiles here, only the final results for the 
electron and spin currents at the depletion layer boundary L, since they respectively determine the charge current in 
the junction and the spin injection through the depletion layer. The boundary conditions and the notation are the 
same as in Sec. ( Ill Al| ). The electron current at L, in analogy with Eq. (|To)) , is J u l — —(D np /L np( ^)F np ^, where the 
modified flux 



F, 



n pC 



Sn L cosh {wp/L n pQ) - 5n p cosh (C,' np w p /2 



5si 



[CnpLnpC 



sinh (wp/Lnpt) 
\ , sinh {C np w p /2) 



p sinh (wp/L n p C )' 



(95) 



If the magnetic drift vanishes, F np Q becomes F np . Since it is J n L which gives the electron contribution to the total 
charge current through the junction, the charge current now explicitly depends on the nonequilibrium spin source 5s p 
and the nonequilibrium spin at the depletion layer boundary, Ssl- These contributions will be important if £ np will 
change on distances smaller than or comparable to L np . Since J„l is sensitive to the sign of C,' np (through the spin 
contribution), the charge current in a magnetic p-n junction could-detect spatial changes in magnetic fields. If the 
junction serves as a solar cell or the base of a junction transistorEj, the nonequilibrium spin 5sl alone will lead to 
charge current, in analogy with the term Sul leading to the usual solar cell current. In fact, both the nonequilibrium 



spin and carrier densities will be normally present when the junction is illuminated by light at x 
of Qnp then either reduces or enhances the solar cell current, depending on the sign of (' np . 
The spin current at L is J s l = —(D np /L np ^)F sp ^, where 



The slope 



5s L cosh (wp/L n p C ) - Ss p cosh ((' np w p /2) 



sinh {w p /L npC ) 



Sn L^ {Cnp L npc) 



— Sn. 



sinh (C n pWp/2) 
'sinh (wp/L npi ) 



(96) 



When neglecting 1/L lp in Eq. (|94|) the equations for electrons and spin become symmetric, so the spin current is 
obtained from the electron current by changing n p to s p and til to sl, and vice versa. Also, in our limit of large L± p , 
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the effective spin diffusion length is L sn — L pn . In this limit the above equation reproduces J s l from Eq. ( |l5| ) . If spin 
relaxation is slow, the spin current in a homogeneously spin split p region does not explicitly depend on the electron 
density. A finite Q' pnl however, couples the electron and the spin densities and the spin current acquires an explicit 
dependence on Sn p and Sn^. Spin injection is modified by magnetic drift too. If the n region remains magnetically 
homogeneous, the injected spin Ssr can be obtained by equating J s l calculated above and J s r from Eq. (pf|). The 
result can be written as 

Ss R = lo Ss n - ( ® np , sn ) sinh (w n /L sn ) F° pC , (97) 



where F® p( - is F sp £ given by Eq. (j9(^) with Shl and Ssl calculated from Eqs. (^2|) and ( |33] ) using San — 0, that is, 
Srii, — ?ioL[exp(V r ) — 1] and 8sl — sql[^p{V) — 1]. Spin injection is modified in several ways. First, there are obvious 
modifications due to changes in the decay lengths, from L np to L np Q. Second, in our limit of \jh\ p = there is no 
explicit contribution of 5n p to 6sr (see Eq. (pq) with 71 =72). Such an explicit dependence appears now because 
of the magnetic drift. Since the factor with Sn p in Eq. (|97j ) changes sign with C,' np , spin injection can be reduced or 
enhanced. Finally, the large bias spin extraction will be affected, since it now depends not only on sql but also on 
?iol- The latter factor again enhances or reduces the large bias spin injection depending on the slope of £ nf) (more 
precisely, on the sign of Cnp s oi)- 

In the opposite limit of fast spin relaxation (which is perhaps more realistic in magnetically doped samples under 
consideration) one can assume for the spin to follow the local carrier density changes: s — a$n. Only the drift-diffusion 
equation for electrons, Eq. (|93]), needs to be solved in this case. To simplify the discussion, we further assume that the 
homogeneous part of the magnetic spin splitting is large, and ao ~ ±1, with a' « 0. The carrier and spin currents 
have the same magnitude, only the sign can differ if «o = — !■ It thus suffices to look at the carrier current. In analogy 
with the previous case, the spin current is determined by F np ^, which now reads 



F npC = Sn L [coth(w p /L npC ) - (aoCnpLnpc) / 2 ] 
cxp (qpC'wp/2) 
p sinh (w p /L np( ) 



(98) 



The spin current and the spin injection (that is, Ssr) are then given as in the previous limit of slow spin relaxation, 
but with F sn £ = aoF np Q. As in the case of slow spin relaxation, here too the I-V curve becomes explicitly dependent 
on magnetic drift. The strength of the magnetic drift is determined by the parameter CnpLnp(, wriue the sign (whether 
it will enhance or reduce the charge current) on the sign of aoQ' np . The solar cell current coming from 5n p depends 
exponentially on . The same applies to spin injection. 



V. SUMMARY AND CONCLUSIONS 



We have studied spin-polarized bipolar transport in magnetic p-n junctions under the general conditions of applied 
bias and externally injected (source) spin. We have introduced a model, by generalizing the successful Shockley model 
of nonmagnetic p-n junctions, to include spin-split bands and nonequilibrium spin. The model is valid only at low 
injection (small biases), although it shows the trends of what to expect at large biases as well. Our theory gives 
the carrier and spin density profiles in the bulk regions (away from the depletion layer), and explicitly formulates 
the boundary conditions for the densities at the depletion layer. In analogy with the original Shockley model we 
employ the condition of (quasi) thermal equilibrium across the depletion layer even when a bias is applied and a 
nonequilibrium spin is injected. However, the spin polarized case requires an additional condition to obtain all the 
relevant input parameters. This condition we formulate in terms of the continuity of the spin current across the 
depletion layer. The obtained boundary conditions allow us to generalize the standard diode formulas to the case 
of spin-polarized magnetic diodes, resulting in a new formulation of the I-V characteristics. Although to explain the 
physics of bipolar spin-polarized transport we use spin polarized electrons only, we also give all the formulas needed 
to calculate the I-V curves for spin polarized holes as well (in the Appendix, where we also show how the equilibrium 
properties of p-n junctions are modified in the presence of spin-split bands). 

We have applied our theory to several cases which we believe are important for spintronics. We demonstrate that 
only nonequilibrium spin can be injected across the depletion layer. Effective spin injection from a magnetic into 
a nonmagnetic region, without a source spin, is not possible at small biases. We show how this claim is relaxed 
at large biases, which build up a nonequilibrium spin in the magnetic majority region, and then inject this spin 
into the nonmagnetic minority region. Similarly, we demonstrate that spin can be extracted at large forward biases 
from the nonmagnetic majority region to the magnetic minority one. We also study spin injection by the minority 
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carriers to the majority region. Physically, this process can be described as spin pumping, since the resulting ac- 
cumulation/amplification of spin in the majority region depends on the spin current of the minority carriers. The 
accumulated spin can be greater than the source spin, which in effect is a spin amplification. A realization of the 
spin- voltaic effect is found at the interface (here the depletion layer) between the minority magnetic region (p) and the 
nonmagnetic but spin-polarized majority region. The spin-voltaic effect is demonstrated by the generation of charge 
current by nonequilibrium spin (at no applied bias). This is also a spin- valve effect, since the direction of the charge 
current can be reversed by reversing an applied magnetic field. The spin induced nonequlibrium charge density is also 
the basis for the spin capacitance of the spin-nalarized junctionsEj as well as for the spin and magnetic field dependent 
charge capacitance of magnetic p-n junctions.Eil Next we have studied (source) spin injection by the biasing electrode 
and shown that this is not a very effective means of spin injection, at least for a simple model considered. Finally, 
we demonstrated that if the neutral regions have nonequilibrium band spin splitting, the resulting magnetic drift can 
significantly affect both the TV characteristics of the junction and the junction spin injection capabilities. 

Our theory is general enough to be applicable to various semiconductor spintronic devices operating under the 
conditions of small injection and nondegenerate carrier statistics. While we have already demonstrated the extensive 
generality of the theory by applying it to a large number of specific model device simulations, we envisage many 
more potential spintronic-jiUnction devices where our models will be useful. Such devices can be, for example, bipolar 
spin junction transistorsEj or spin thyristors, with great technological potentials, and where charge currents (and 
their amplification) can be controlled not only by bias, but also by nonequilibrium spin and magnetic field. However, 
to apply the theory to realistic device structures, many physical aspects of the model will need to be modified. In 
many cases the spin states of the carriers are not simple spin doublets, but rather multiplets, as a result of the 
spin-orbit coupling. In addition, the electron-hole recombination is, in general, spin selective, so if both electrons 
and holes are spin polarized, more realistic models for the recombination need to be introduced. Furthermore, 
carrier recombination and spin relaxation depend on the carrier density, an effect which may be found important if 
ferromagnetic semiconductors are employed. Other possible additions to the model may include a realistic treatment 
of spin relaxation (and carrier recombination) in the depletion layer and finite spin relaxation at the contact electrodes. 
Structural modifications may include inhomogencous magnetic doping (or inhomogeneous magnetic fields) also in the 
bulk regions, and schemes based on two or three dimensional spin bipolar transport. Since, at the moment, there is a 
lack of experimental understanding of bipolar spin transport, theoretical modeling (both analytical as presented here 
or numerical, which is of greater applicability, as reported in Refs. |52] , |53| , |54| ) is particularly important. We believe 
that although quantitative aspects of spin-polarized bipolar transport may be seriously modified, our theory captures 
the essential physics and the predicted phenomena are robust enough to be present in more realistic situations. 



Acknowledgments 



We would like to thank Prof. E. I. Rashba for useful discussions. This work was supported by DARPA, the US 
ONR, and the NSF. 



APPENDIX A: EQUILIBRIUM PROPERTIES OF MAGNETIC P-N JUNCTIONS 

To study equilibrium properties of magnetic p-n junctions we consider both electrons and holes spin polarized. 
Denote the electron and hole spin densities as s n and s p , and reserve the second subscript (if needed) to denote the 
region. Symbol denotes the equilibrium values. As in the main text, the energies (potentials) are given in the units 
of UbT {ksT/q). Further, denote as 2£ ra and 2( p the spin band splittings of the conduction and valence bands. We 
adopt the convention that £„ (£ p ) is positive when the spin up electrons (holes) have a lower energy than those in the 
spin down states. Both electrons and holes are assumed in thermal equilibrium, obeying nondegenerate Boltzmann 
statistics. 

The equilibrium carrier densities obey the law of mass action, now reading 

n p Q = n\ cosh(C„) cosh(C p ). (Al) 
As a result, the minority carrier densities (electrons in the p region and holes in the n region) are 

o 
71- 

n 0p = -^-cosh(C„ p )cosh(C )9 p), (A2) 

9 

71- 

Pon = -T7-cosh(C P „)cosh(C„„). (A3) 
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Similarly, the corresponding equilibrium spin densities are 

2 
TV? 

So, np = -^sinh(C„p)cosh(C pp ), (A4) 
n 2 

so, P n = -77- sinh(Cpn) cosh((„„), (A5) 

so that the equilibrium spin polarizations of electrons and holes are ao n = tanh(£ n ) and uq p = tanh(C p ). The built-in 
voltage, which is the electrostatic potential drop across the depletion layer depends on the band splittings (thus on 
the equilibrium spin polarizations): 

V b = In f -^r^ J - ln [cosh(Cpp) cosh(C„„)] ■ ( A6 ) 

The built-in voltage is slightly reduced by the spin splitting. Note that only the band splittings of the majority 
carriers affect the built-in field. The reason is that the chemical potentials in the bulk regions (considered separately) 
depend only on £„„ and ( pp : 

Mon = Mi + m f — J - lncosh(C„„), (A7) 

Hop = P-i - In + hicosh(Cpp). (A8) 

Here fj,i is the chemical potential for the intrinsic (undoped) and unpolarized case. In making a junction, the built-in 
field arises upon equilibrating the two chemical potentials: Vb = fio„ — fj,o p . The band splitting does not affect the 
nondegeneracy of the carrier statistics, since the distance between the chemical potential and the lower conduction 
(upper valence) spin band does not change with £ at large C- 



APPENDIX B: SPIN POLARIZED HOLES 



Spin polarization of holes can be treated separately from that of electrons, since, in our model, electron and 
hole transport are independent (only minority diffusion is considered), and the electron- hole recombination is spin 
independent (in our simplified picture electrons of one spin can recombine with holes of either spin). Inclusion of 
spin polarization of holes into our theory then amounts to simple notation exchange, p with n and L with R. For 
completeness, we present all the important formulas which arc needed to obtain the charge current contribution by 
spin polarized holes. Since this is a separate section from the main text, we adopt the same notation for the hole 
spin as we had before for electrons, without using more elaborate set of indexes. The hole spin density (only in this 
section) is s and the hole spin polarization is a. All the other symbols retain their original meaning. 

In analogy with Eq. (Ell), the hole charge current is 



jp = jo P + hv + hv, 



(Bl) 



where 



1) 



JOp — 3gp ( e 



V r a 0R - ®0L 

Jip = J 9P e da L - 2 — . 

1 a 0L 



1 



J2p — ]gp 



5p„ 



The hole generation current is 



Jgp 



qD pr . 



cosh(w„/L p „) pon 



POR COth {Wn/L pn ) . 



The (now majority) hole spin, in analogy with Eq. ( p6| ) can be expressed as 

Ss L = j Ssp + 7i<5s„ + j 2 a°5p. 



13S 0R (e v - 1) , 



(B2) 
(B3) 

(B4) 



(B5) 



(B6) 
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where the geometric/transport factors are 



70 = 1/ cosh(w p /L sp ), (B7) 



D pp J \L sn J smh(w n /L sn )' 
D pn \ ( L sp \ tanh(w p /L S p) 



72 = ^ ;;?7 it y (B9) 

\DppJ \L pn J smh(w n /Lp n ) 

73 = 72COsh(u)„ /L pn ). (BIO) 
Finally, the (now minority) hole density and spin in the n side of the depletion layer are 

v f, , .otoR-aoL 



'0L 



PR = PoRe [l + 6a L — g — ' ( B11 ) 

1 - at 



v i Sa L 1 - aoRao L \ 

SR = SQR6 1 H 5 ' ( B12 ) 



aoR 1 



Physical consequences of the spin polarization of holes in bipolar transport are in complete analogy with the physics 
discussed in t he main text where only spin-polarized electrons are considered. In particular, the hole charge current j p 
from Eq. flBl| ) needs to be substituted to the total charge current formula, Eq. ([!?]). In many cases, however, one can 
realistically treat only one carrier type as spin polarized. If, for example, holes have a very short spin lifetime (or small 
diffusivity) , their spin polarization (even their contribution per se) does not need to be considered. The e xcepti onal 



cases are the large bias spin-polarized transport and magnetic drift in the bulk regions, treated in Sees. IV E and 



IV F, respectively, in which electron and hole transport can be strongly coupled. 



APPENDIX C: MAJORITY ELECTRON DRIFT AND DIFFUSION 



The spin profile in the n region is affected by the electr ic field and charge neutrality, as described by Eq. (|7£ 
Assuming the same boundary conditions for spin as in Sec. Ill A 2| , that is, 5s(d n ) = 5s r and 5s(w n ) = Ss n , and the 
boundary conditions for holes 8p(d n ) — 6pr, and Sp(w n ) = (Ohmic contact), the solution to Eq. (|79| ) can be written 
in the form analogous to Eq. (18): 



Ss 



-E(x-d n )/2 



(5s R coshrjsE + F sE sinh?7 s£ ) + a 0n 5p. 



(CI) 



We now describe the new notation. The effective deviations from the equilibrium of spin and hole densities are 

Ss R = Ss R - Aa Qn SpR [1 + X E coth (w n /L pn )] , (C2) 
Sp = A[5p + X E {L pn /D pn )Jp], (C3) 

where Xe = E/ n 2 L pn measures the strength of the electric field for drifting spin, k = \/(V-kpn — V-^sn) (the- singular 
case of k = is excluded from the solution), and 1/A = —(L sn K) 2 (l — X E ). The normalized flux is 



F S E = 



-5sRCOsh(w n /L sE ) + Ss n exp (Ew n /2) 
sink (w n /L sE ) 



(C4) 



introducing a length scale L s e for electric spin drift: 1/L s e = \f[E 2 / 4 + 1/L 2 n ), which is also used to define 
r] S E = (x — d n )/L S E- This length scale was already introduced in Refs. |55|,p6|,p4]. The new effective spin density 



5§ n = Ss n - 



Aa 0n 6p R XE 
sinh (w n /L pn )' 



(C5) 



Finally, the hole density profile 5p is obtained by solving independently for hole diffusion, Eq. ( p0| ) . For completeness 
we show the result: 



where r\ pn = {x — d n )/L pn and 



Sp = Spr cosh r\ vn + F pn sinh rj pn 



-Sp R coth (Wn/L pn ) 



(C6) 
(C7) 
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The hole current is then J p = —D pn Sp'. 

The importance of electric drift for the majority electron spin transport is in aiding the spin injection through the 
depletion layer, from the majority magnetic to the minority nonmagnetic region. To see how spin can be injected 
through the depletion layer we need to know the s pin current at the depletion layer boundary. The spin current profile 
is J s = —D nn (sE + s'), where s is given by Eq. (CI). The spin current at x — d n is 



J, 



sR 



1 



t(so« + -Ss R )E 



Erin \ 

Am 5s R cosh (w n /L sE ) - Ss n exp (Ew n /2) 



D 



sinh(w n / L s e) 



+ — a^ n 8p R A coth ( — — 



L 



w„ 



(C8) 



^^-a 0n Sp R A coth 

^sE \^sE 



where we neglected terms of order \ E E nn bp R j L pn . For most practical cases in magnetic p-n junctions the electric 
field at low injection can be neglected, so that EL pni EL sn <C 1. Then L s e ~ L sn , k rj 1/L pn , and A w — (L pn /L sn ) 2 . 
Since E is of order 5p R /L pn , the contribution to the spin current (and thus to spin injection) from the hole density 
is negligible, since in the considered limit A <C 1. If Ss n is greater than, say, 10~ 3 x N^, then also the contribution 
from the electric drift can be neglected (not limited to the above limit), verifying our theory in the main text. If, 
however, the source spin is small, and there is appreciable forward bias, the electric drift has to be taken into account 
for describing spin injection across the depletion layer. This is done in Sec. (IV E). 
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